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^ We describe an approximate quantum mechanical picture of hadrons 

K-5 in Minkowski space in the context of a renormalization group procedure 

for effective particles (RGPEP) in a light-front Hamiltonian formulation 
of QCD. The picture suggests that harmonic oscillator potentials for con- 
stituent quarks in lightest mesons and baryons may result from the gluon 
' ^ ' condensation inside hadrons, rather than from an omnipresent gluon con- 

densate in vacuum. The resulting boost-invariant constituent dynamics 
2^ at the renormalization group momentum scales comparable with Aqcd, 

is identified using gauge symmetry and a crude mean-field approximation 
for gluons. Besides constituent quark models, the resulting picture also 
I— I resembles models based on AdS/QCD ideas. However, our hypothetical 

picture significantly differs from the models by the available option for a 
^ systematic analysis in QCD, in which the new picture may be treated as a 

candidate for a first approximation. This option is outlined by embedding 
our presentation of the crude and simple hadron picture in the context of 
I RGPEP and a brief outlook on hadron phenomenology. Several appendices 

\^ describe elements of the formalism required for actual calculations in QCD, 

including an extension of RGPEP beyond perturbation theory. 

o 
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1. Introduction 

Contrast between the simplicity of constituent quark model (CQM) clas- 
sification scheme for hadrons in particle data tables [1] and the complexity 
of QCD, exists since quarks were proposed to explain hadrons [2] and QCD 
was proposed as a theory of strong interactions [31 IH [5] . The contrast can be 
illustrated by the simplicity of attempts to build relativistic quark models, 
in which one introduces harmonic oscillator potentials [6], in comparison 
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with the complexity of QCD sum rules j7] , in which even the state without 
any hadron, i.e., vacuum, is meant to contain a complex structure. The 
complex structure is reflected in the concept of condensates. 

The vacuum structure in quantum field theory (QFT) eludes theorists [5] 
despite that they have a complete command of perturbative calculations. 
For example, in the case of QCD, theorists can use asymptotic freedom 
[HI Ho] to calculate quark and gluon scattering amplitudes but they also have 
to use a parton model [TT] to relate their calculations to observables. The 
point is that the mechanism of binding of partons continues to be a mystery. 
The mystery is understood in principle as related to the complex features 
of the theory that lie beyond reach of perturbative approaches. Regarding 
these complex features, palpable progress in calculating hadronic masses is 
achieved using lattice formulation of QCD [121 US]- However, so far it is 
not clear how to use QCD for producing quark and gluon wave functions of 
hadrons in the Minkowski space and thus provide a picture of hadrons of 
comparable precision to the picture of atoms achieved in QED. 

In this situation, it is interesting to note that there exists a formulation of 
QCD which uses Dirac's light-front (LF) form of Hamiltonian dynamics |14[ 
[T5] and does not introduce any complex vacuum structure [16j. How could 
then the LF formulation of QCD produce the effects that in other approaches 
are associated with the concept of a complex vacuum state, including effects 
associated with spontaneous symmetry breaking or the gluon condensate? 
Ref. [16] points out that an effective Hamiltonian expected to result from 
a renormalization group analysis of the canonical LF Hamiltonian of QCD 
with all counterterms that are necessary for stabilizing the analysis, may 
contain new terms that provide such effects. 

Here it is argued that the renormalization group procedure for effective 
particles (RGPEP, see Section [2]) provides a new way of thinking about 
light quarks in LF QCD. In this new way, the harmonic oscillator poten- 
tials, previously associated in Ref. [L7j with a gluon condensate in vacuum, 
can be interpreted as coming from the gluon content of a hadron rather 
than a vacuum. Previous reasoning in Ref. [17] was developed using the 
instant form pj:] of dynamics, i.e., the standard Hamiltonian evolution in 
time. Now, the possibility of considering harmonic oscillator potentials for 
light quarks at distances comparable with the characteristic size of hadrons, 
meant here to be equivalent to the distances on the order of 1/Aqcd in the 
RGPEP scheme, is pointed out using a decomposition of hadronic states 
into the LF Fock space components, i.e., states obtained from empty vac- 
uum using creation operators defined in the LF quantization of fields rather 
than the standard canonical quantization of the instant form. These Fock 
components contain various numbers of virtual effective particles whose in- 
teractions depend on the renormalization group scale in the RGPEP scheme. 
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The scale is denoted by A. One can think about A as having dimension of 
momentum. Low-energy features of hght states in QCD made of Ught quarks 
are expected describable using an effective Hamiltonian that corresponds to 
A comparable with Aqcd- Readers interested in details of RGPEP that 
are already known in the case of heavy quark dynamics, can consult Ref. 
|18| , where an approximate QCD calculation of masses of heavy quarkonia 
is described. Here only the case of light quarks is discussed. 

According to Ref. [T^, if gluon condensate effects are coming from the 
gluon content of hadrons rather than vacuuirl^ one can avoid the problem 
of an excessively large vacuum energy density in cosmology. However, the 
new reasoning presented here does not deal with cosmological issues. This 
article concerns only the particular mechanism that was described in Ref. 
|17j . In the mechanism in Ref. [T7j, the gluon condensate meant to be in 
vacuum provided harmonic oscillator potentials for constituent quarks in 
light mesons and baryons. The potentials agreed with known phenomeno- 
logical CQMs j211 [22l [23] when the expectation value of the gluon field 
strength squared was equated to the vacuum gluon condensate value used 
in the QCD sum rules [[7]. It is argued below that it was not necessary 
to assume, as it was done in Ref. [T^, that the gluonic expectation value 
came from the empty vacuum. The point of view developed here is that 
the gluon expectation value might be coming from the content of a hadron 
itself, and the relevant content can in principle be identified in LF QCD us- 
ing RGPEP. Potentially broad implications of such change in interpretation 
are not discussed here. The hypothesis put forward in this article requires 
a considerable amount of work to verify and it would be premature to draw 
broad conclusions before relevant RGPEP calculations are completed. 

The calculations described here arrive at similar harmonic potentials to 
those found in Ref. pT] but in a different way and with a significantly dif- 
ferent interpretation. Namely, instead of starting from the non-relativistic 
(NR) Schrodinger equation and gluon vacuum condensate, one starts from 
the LF QCD Hamiltonian eigenvalue equation for mesons and baryons. 
Canonical quark and gluon field operators are formally transformed using 
RGPEP to effective operators at the running cutoff scales that are compa- 
rable with Aqcd- Using color gauge invariance principle, together with a 
mean-field approximation that still needs justification and is not verified yet 
by a rigorous ab initio calculation, one arrives at the eigenvalue equations in 
which expectation values of effective gluon fields are evaluated in a gluonic 
component of hadronic states rather than in the vacuum. The gluonic com- 
ponent is postulated to be universal. Then, using new relative momentum 
variables, one obtains LF wave functions that describe the relative motion 



^ The quark condensate has also been suggested to originate in hadronic content instead 
of vacuum [20] . 
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of hadronic effective constituents in a boost-invariant way. In summary, the 
reasoning described here imphes a relativistic picture of hadrons in which 
the same harmonic oscihator potentials that agree with CQMs phenomenol- 
ogy may result from the gluon content of hadrons in LF QCD rather than 
from the vacuum per se. 

It should be clarified that harmonic potentials among effective con- 
stituents in lightest hadrons apply in a situation where color charges are 
close to each other. If they were separated by a distance considerably larger 
than the size of a light hadron, their invariant mass would increase con- 
siderably due to the potential. Interactions other than the potential term 
alone would be activated, able to create additional particles and changing 
the dynamics, perhaps including a string of effective gluons. 

The article is organized in the following way. Section [2] describes the 
concept of effective particle in the context of RGPEP and explains why 
the vertex form factors (denoted by fx) eliminate changes of the number of 
massive virtual effective particles. This feature is required for thinking that 
a hadron can be represented by a convergent expansion in the basis of the 
Fock space. The basis that counts is not built in terms of canonical quark 
and gluon operators but in terms of the effective ones, corresponding to 
relatively small scale parameter A. The convergence is not proven, but it is 
deemed not excluded given the exponential falloff of the RGPEP form fac- 
tors fx as functions of constituent invariant masses. Section [3] outlines the 
representation of states of light hadrons in the effective particle basis in the 
Fock space. Section |4] discusses a perturbative expansion for the unitary 
transformation Ux that connects current, or canonical quarks and gluons 
with their CQM counterparts, and for the transformations WxiX2 that con- 
nect effective particles that correspond to different values of the scale A. The 
eigenvalue problem for light hadrons is considered in Section [5j This Section 
describes our derivation of LF invariant mass operators for quarks including 
their minimal coupling to the gluons condensed in a hadron. Next Section 
[6] shows how the RGPEP can be applied in phenomenology of hadrons, 
with emphasis on form factors, structure functions, and connection with 
AdS/QCD ideas. Section [T] concludes the paper. Finally, four appendices 
are provided in order to support the claim that RGPEP can be used to verify 
our reinterpretation of the gluon condensate in QCD. Appendix \K\ outlines 
how CQMs can be viewed as limited representatives of the same universality 
class that QCD belongs to. Appendix [B] describes details concerning per- 
turbative calculations of RGPEP transformations. Presentation of RGPEP 
beyond perturbation theory is given in Appendix [Cj The last Appendix [D] 
offers a visualization of the RGPEP scale dependence of hadronic structure. 
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2. Effective particles and vertex form factors 

Reasoning described in tlie next sections requires elements of RGPEP 
[21] This section focuses on the concept of effective particles and the role 
of vertex form factors that RGPEP introduces in the interactions of effective 
particles, starting from QCD. 

In QCD, RGPEP begins with the canonical LF Hamiltonian in which 
the dynamically independent quark and gluon quantum fields are expanded 
into their Fourier components. These components are the operators that 
create or annihilate single bare quarks or gluons of definite momentum. For 
brevity, creation and annihilation operators will be commonly called particle 
operators, in order distinguish them where necessary from quark and gluon 
field operators. 

The canonical LF Hamiltonian contains singular interaction terms. It 
must be regulated to avoid infinities. This is done by introducing regular- 
ization factors in interaction terms. These factors are constructed to limit 
relative motion of particles that participate in interactions. For example, let 
the total momentum of all particles in an interaction term have components 

and P"*". The notation means: = P^ + P^, z-axis is the direction 
distinguished in the definitior|^of the LF, _L denotes transverse components, 
i.e., the components x and y that are transverse to the z-axis. Let a selected 
particle have momentum components p"*" = xP^ and p'^ = x P-^+k-^. Then, 
the regulating factor for this particle that is sufficient for taming logarith- 
mic divergences in QCD can be of the form x'^exp {—k-^ / IS?) [25], where 
(5 — )• and A — )• oo are the small- j; and ultraviolet regularization param- 
eters, respectively. Initially, one also introduces an absolute infinitesimal 
lower bound e"*" < p"*" for all particle operators in order to eliminate the LF 
zero modes. This step amounts, however, to saying that all particles must 
have positive momentum component p"*". This requirement eliminates com- 
plex vacuum. Subsequently, one demands that for every particle in every 
interaction term x > e, with an infinitesimal number e. The number e is a 
fraction of momentum. It is not related to e"*". Then, when one works in the 
limit ex^^ — t- point-wise in x when e is sent to 0, the factor x^ takes over 
the regulatory function from e. As a result, the variables x range from to 
1 and there are regulating factors x^ and (1 — x)^ at the end points in all 
interaction vertices. The variable x = p'^ / P^ for a particle of momentum p 

^ Appendix [a] describes RGPEP in the context of asking how QCD can be transformed 
into an effective theory that resembles CQMs as approximate representatives of the 
same universality class. 

LF is a hyper-plane in space-time that is swept by a wave front of a plane wave 
of light. In standard notation, the frame of reference of the inertial observer who 
develops a quantum theory is set up so that the plane wave moves against the z-sxis 
in the chosen frame. 
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in a hadron of momentum Ph coincides with the parton model longitudinal 
momentum fraction in a hadron in the infinite momentum frame (IMF). 
Variable /c"*" is equal to the transverse momentum of a parton in the IMF 
when the parent hadron of the parton moves precisely along z-axis. Note, 
however, that LF QCD is formulated in terms of the same variables x and 
fc-*- no matter how hadrons move, i.e., the LF Hamiltonian theory is explic- 
itly invariant with respect to the 7 Poincare transformations that preserve 
the LF. This means that the theory of a hadron structure in a rest frame of 
the hadron has the same form as in the IMF. 

When the regularization is being removed, 5 — )• and A — )• oo, RGPEP 
establishes the required ultraviolet counterterms [2l[j. The small-x regu- 
larization drops out from dynamics of colorless states because they cannot 
produce long-distance interactions along the LF. 

Let the regularized LF Hamiltonian of QCD with counterterms be de- 
noted by H and let b denote bare particle operators. RGPEP introduces 
effective particles of scale A through a unitary transformation, 

bx = UxbUl (1) 

The corresponding Hamiltonian operator, 

Hxibx) = H{b), (2) 

is a combination of products of operators bx with coefficients cx that are 
different from coefficients c of corresponding products of operators b in the 
canonical Hamiltonian with counterterms, H{b). Since A is related to an up- 
per limit on momentum transfers in interactions, the operators b correspond 
to A = oo. For the infinite A, we have b^o = b and H{b) = -ffoo(&oo)0 

RGPEP provides differential equations that produce expressions for the 
coefficients cx in Hx- The operator Ux is calculable in RGPEP order-by- 
order in an effective coupling constant [24J in the form of normal-ordered 
products of particle operators bx ■ Appendix [B] illustrates how it is done in 
lowest orders. Appendix [C] shows how non-perturbative calculations can be 
attempted. Appendices [A| |B} and [C} provide a formal background for the 
entire discussion that follows. 

The key feature of Hx{bx) that emerges from RGPEP is that all coeffi- 
cients Cx contain a form factor fx. The form factor prevents changes of the 
total invariant mass of effective particles undergoing interaction by more 
than A. The RGPEP parameter A plays thus the role of a momentum-space 
width of vertex form factors in all interaction terms. In fact, RGPEP is 

* Where it is unlikely to lead to a confusion, Hx{b\) is abbreviated to Hx- Later, 
also operators such as H\^{bx2) occur, meaning an operator that is a combination of 
products of particle operators bx^ with coefficients cx^ instead of c^j ■ 
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designed to work this way. When A is smah, the form factors suppress all 
interactions among massive particles besides the terms one can call poten- 
tials, i.e., the interaction terms that do not change the number of massive 
particles (see below). In other words, the RGPEP transformation Ux is de- 
signed to identify the dynamical relationship between nearly massless cur- 
rent quarks and their field-theoretic interactions at formally infinite A with 
massive constituent quarks and their interactions through potentials at A 
so small that potential models may apply as an approximation to solutions 
of the whole theory for states of smalles masses. 

If one so desires, transformation U\ can be kinematically complemented 
with the Melosh transformation [26j , in order to express the spinors typically 
used for description of constituent quarks in the constituents fram^ (CRF) 
with the spinors one can conveniently use in the IMF0 

It follows from its definition that the form factor fx prevents the num- 
ber of effective particles from changing if their masses m\ exceed A. For 
example, suppose that a virtual particle of mass mx emits another virtual 
particle with mass mx- The change of invariant mass of particles in the 
interaction is at least mx- The form factor fx is designed to quickly tend 
to for invariant-mass changes greater than A. So, if mx exceeds A, the 
form factor eliminates the emission. In particular, if effective gluons are 
assumed to have effective masses nix, the Hamiltonian Hx will not include 
interactions that can change the number of effective gluons when nix > 
A sizeable mass of effective gluons in Hx is what one may expect to happen 
at small values of A in QCD because such effective gluon mass is a candi- 
date for explaining the absence of small spacing in the spectrum of hadronic 
masses. There should be small spacing in the presence of massless gluons, as 
it happens in atoms described by QED with massless photons. In contrast 



^ The CRF frame is an inertial reference frame in which the total momentum of con- 
stituents treated as free particles of definite masses has only time component different 
from zero. The concept of CRF in LF dynamics is different from a similar concept in 
the instant form of dynamics, because conservation of in interactions makes the 
CRF differ from the center-of-mass system (CMS) of a hadron in the LF dynamics, 
while the CRF and CMS are the same in the instant form of dynamics, where is 
conserved by interactions. P^ is preserved in interactions in both forms of dynamics 
equally. 

^ Ref. [4] is of interest here because the last sentence in it suggests one might perhaps 
use some collective coordinates as a satisfactory method for truncating QCD, instead 
of "the brute-force lattice gauge theory approximation!" While Ref. @] is quite new 
to the author at the time of writing this article, it should be noted that RGPEP 
can be viewed as an attempt of this type. Namely, the effective quarks and gluons, 
as constituents of size ~ ^qcd' describe collective modes in dynamics of many 
small quarks, anti-quarks, and gluons that nevertheless may be individually active in 
virtual processes characterized by large invariant-mass changes, allowed in a single 
interaction only when A ^ Aqcd- 
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to atoms, hadrons do not exhibit such small spacing. 

Note that a smooth form factor fx must allow for some small, transi- 
tional range of values of A < mx. In the transitional range, interactions 
that change the number of effective gluons gradually disappear when A is 
lowered below mx- In the reversed RGPEP evolution in A from small to 
large values, gluons will gradually appear in the dynamics as A increases 
and mx decreases so that at some point A becomes greater than mx- 



3. RGPEP representation of states of light hadrons 

If QCD is to explain precisely the success of classification of light hadrons 
in terms of constituent quarks with masses on the order of 1 /3 of a nucleon 
mass, a clearly defined procedure must explain how the masses m„ and of 
the lightest quarks increase from their standard model (SM) values of order 
5 MeV in a local gauge theory to their common constituent value of order 
Aqcd in a corresponding effective theoryj^ It is not known yet if RGPEP 
leads to such increase of effective masses of light quarks in QCD when the 
parameter A is lowered toward Aqcd, or even below Aqcd- We assume 
here that RGPEP does lead to such result]^ This assumption cannot be 
rigorously verified yet because the domain of A ~ Aqcd cannot be reached 
using low orders of perturbation theory and non-perturbative solutions to 
RGPEP equations are still not known in QCD. 

When A is large in comparison to Aqcd, so that the effective coupling 
gx is small and perturbation theory applies, the eigenvalue problem in- 
volves many Fock sectors with many quarks of small Lagrangian masses. In 
these circumstances, it is hard to analyze the eigenvalue equations precisely. 
When A is lowered, the number of necessary Fock sectors is expected to de- 
crease because of the vertex form factors fx but the coupling constant gx 
increases and it is hard to evaluate Hx precisely. A way out of this situation 
is to calculate Hx in terms of successive approximations]^ 

A workable RGPEP method for building successive approximations has 
been outlined including results for hadron masses and wave functions in the 
case of heavy quarkonia [18j . where the number of heavy quarks is nearly 
constant, equal 2. In that case, one can take advantage of the NR approx- 
imation for relative motion of quarks, based on smallness of the effective 

^ How RGPEP is able to relate a local canonical theory to a non-local effective theory 
is explained in [27J. 

* Quark masses increase as momentum scale is lowered in perturbation theory de- 
veloped in terms of the Feynman diagrams 28 . Similar results are obtained from 
Dyson-Schwinger equations [25]. LF constituent models can incorporate running 
masses [50] , 

^ This includes successive approximations for solutions of the non-perturbative RGPEP 
equations in Appendix [Cj 
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coupling constant gx when A is comparable with heavy quark masses that 
are much larger than Aqqd- 

However, no counterpart for such approximation scheme is available yet 
for light quarks in LF QCD using RGPEP. The reason, stated already above, 
is that, in distinction from the case of heavy quarks, the light quark masses 
are so much smaller than Aqcd that many light quarks can a priori be 
created or annihilated in the interactions that limit virtual energy changes 
by about Aqcd- Therefore, one needs to approximate reasonably well a 
great deal of relativistic dynamics of virtual particles that are constantly 
created and annihilated by strong interactions. It is suggested below that 
a reinterpretation of the gluon condensate as a part of a hadron provides a 
guideline for moving in that direction. 

It is a conceptual jump that requires a connection between constituent 



quarks and structure functions (see Section 6.2), but one can observe here 
that parton distribution functions at small Q"' are typically very stable as 
functions of and quarks always carry only about half of the nucleon 
momentum [31^ [H^ 1551 135) . The gluon component of hadrons, which 
plays the role of a gluon condensate in the picture described below, may only 
be responsible for a part of the other half of hadron momentum. Namely, 
the constituent quarks also contain gluons. The issue of how heavy is the 
part of a hadron that is made of effective quarks at A < Aqcd and how 
heavy is the corresponding gluon condensate component, is discussed in 
Section [6] and Appendix [D} Firm conclusions require RGPEP calculations 
that have not been done yet. 

In the picture described below, the masses of effective quarks and gluons 
at some scale A < Aqcd are assumed greater than A itself. All interactions 
that change the number of the effective particles are absent and the only 
interactions left in Hx with A near and below Aqcd are potentials. These 
potentials are expected to match potentials used in the CQMs. Our rein- 
terpretation of the gluon condensate is associated with a definite result for 
the potentials that act among the constituent quarks. In summary, the 
LF Hamiltonian of QCD with some A = Ac ^ Aqcd is assumed to describe 
quarks of mass rric > Ac. These quarks interact only through potentials, i.e., 
interactions that do not change the number of constituents. The subscript 
c refers to the word constituent. 

An eigenstate \ip) of the effective Hamiltonian Hx^{bx^) that represents a 
hadron of momentum P/j, is a superposition of effective-particle basis states 
(n = 2 for mesons and n = 3 for baryons) 



|l...n)A. = H^IaJO), (3) 

i=l 
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1^) = ^V'Ae(l,-,n)|l...n)A,. (4) 

l...n 

The arguments 1 to n of the wave function provide a shorthand notation for 
momenta, spins, flavors, and colors of the corresponding particles. For ex- 
ample, 2 as an argument of a wave function stands for the thrcc-momcntum, 
spin projection on z-axis, flavor, and color of the particle number 2. Sum- 
mation over numbers 1 to n is a shorthand notation for summation over the 
quantum numbers of the corresponding particles including integration over 
their momenta. In this abbreviated notation, the eigenvalue equation for a 
mass and a wave function of a hadron with momentum and built 
from n constituents of scale Ac, takes the form 

,^{l...n\{Pt^H,^-Pft'm = M',^{l...nm. (5) 

This eigenvalue equation with n = 2 or n = 3 is meant to correspond to the 
CQM picture of hadrons as built from 2 or 3 constituent quarks. 

Note that the vacuum state |0) is not changed when A changes. This is a 
unique feature of LF dynamics. In the standard, instant form of dynamics, 
a change in particle operators must be accompanied with a change of the 
corresponding vacuum state. 

One can write the same eigenvalue equation for the same states using 
other values of A than Ac. The eigenstate is the same for all values of 
A but the basis states and corresponding wave functions depend on A. For 
A S> ^QCDi the state |V') may contain a giant number of multi-particle Fock 
components built by acting with creation operators h^^ on the vacuum, while 
when A = Ac, the entire state contains only 2 or 3 constituent quarks. 

Using RGPEP, one can express constituent quarks at scale Ac in terms 
of effective quarks and gluons corresponding to A > Ac, 

6a, = WhxW'^ , (6) 
W = Ux^ul (7) 

Momenta, spins, and izospins of the effective particles are the same, irre- 
spective of the change in parameter A. Since W^|0) = PF^IO) = |0), one 
has 

n 

|1-")a. = Wj{h\^\Q), (8) 
and the same eigenstate can be written as 



IV') = ^VAe(l,-,n)T^|l...n)A. 

l...n 



(9) 
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The result of action of W on the basis state with n constituents correspond- 
ing to scale A, with A > Ac, is a coherent slew of Fock sectors with various 
particle numbers, with the number of quarks grater or equal to the minimal 
constituent number n for a hadron. The momentum space wave functions 
of the resulting Fock components are determined by the bound-state eigen- 
value condition of H\. 

Although a precise structure of Hx is beyond insight of perturbation 
theory for A ~ ^qcd and W at such scales also contains non-perturbative 
dynamics, some generic features of the Fock wave functions in Eq. ^ can 
be assessed on the basis of generic features of W implied by properties of 
operators and Ux that are visible already in perturbation theory [36j. 
Operator W conserves momentum. W can replace a single effective particle 
with a bunch of other particles which carry together the same momentum 
as the replaced particle. W can also annihilate a whole set of particles and 
create a new set with the same momentum. 

Hence, the colorless meson (M) and baryon (B) states \tjj) that corre- 
spond to the CQM picture (color factors are written explicitly). 



l^>M = ^Va.(1,2)^6^1^4J0), (10) 
12 ^3 



^abc 



)B = T.^X.{W)—bli/lb^,JO), (11) 



123 



are equal to, respectively, 



sab 

5;VA.(l,2)^T^6^14i|0), (12) 



^abc 



)B = J] Va.(1,2,3) ^ Wblib'^b^.lO) . (13) 

123 ^6 



The representation of meson and baryon states in Eqs. (10) and (11) in 
terms of only 2 or 3 effective constituents at scale Ac, is thus transformed 
into the representation of the same states in Eqs. ( 12 ) and ( |13[ ) that are built 
from effective particles of scale A. The resulting Fock space wave functions 
for effective particles of size A~"^ involve functions ipx^il, 2) and ipx^il, 2, 3) 
through convolutions that emerge from the sums in Eqs. (12) and (13) and 
the structure of W. This is how the complex Fock-space structure of QCD 
hadrons is supposed to contain information about a simple CQM picture. 



12 
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4. Structure of W 

In order to provide an illustration of how W may act on quark and gluon 
states, this Section explains the structure of W obtained just in first-order 
perturbation theory in the effective coupling constant g\. The perturbative 
expression is obtained starting from Eqs. ([T| and (A.l), which produce 
together 

= [nx,r] , (14) 

where T = l^ljx^x and Ux = Ux{b^). This means that the functional J-x 



in Eq. (A.l) is set to have the form 

FxlUx] = [nx,T] , (15) 



and the key to RGPEP is the dependence of T on Tlx [21]. This dependence 
is designed keeping in mind that the resulting Hamiltonian Hx{bx) should 
contain the form factor fx in each and every interaction term. Suitable 
notation for this condition is provided by writing Hx = fxGx, where Hx 
contains the form factor fx in interaction vertices while Gx does not. If an 
operator Gx contains a product of creation and annihilation operators bx 
with a coefficient cx, the operator Hx = fxGx contains exactly the same 
product with coefficient fxcx-, where fx depends on the difference between 
the invariant mass squared of particles annihilated by annihilation opera- 
tors in the product and the invariant mass squared of particles created by 
creation operators in the product. Both invariant masses are calculated 
using the particle kinematical momentum variables and eigenvalues of part 
Hqx = Gqx of the Hamiltonian Hx- 

Thus, using Hx = Hx{b) and omitting A, one writes 

n = fg = go + fQi. (le) 

Differentiating % with respect to A one arrives at the equation that defines T 
in such a way that it must vanish when interactions vanish (this guarantees 
that T is expandable in powers of the coupling constant) and that its matrix 
elements between eigenstates of Qqx vanish when the differences between the 
corresponding eigenvalues vanish 



[TM = [(i-mr. (17) 

Prime denotes differentiation with respect to A. The solution is denoted by 

r = {[{l-f)9l^]g^ , (18) 
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where the curly braces with subscript Qq indicate the energy difference (dif- 
ference of eigenvalues of Qq) in denominator that results from the commu- 
tator of T with Qq on the left-hand side in Eq. (|17|). Then lA satisfies 



w = z^{[(i-mr}c„ , (19) 



with condition Uoo = 1. Eq. ( 19 ) can be expanded in powers of the coupling 
constant. 

^.1. Perturbative expansion for Ux 

Let us assume here that Qq does not depend on the coupling constant 
to all orders of perturbation theory. This means that perturbative self- 
interaction terms are included in Qj and Qq is actually independent of A 



10 



While this condition may seem quite restrictive, note that it allows for 
appearance of mass terms that are proportional to positive powers of Aqcd, 
since Aqq d vanishes to all orders of perturbation theory. Writing expansions 
in the bare coupling constant, 

Ux = 1 + guxi + 9'^ux2 + (20) 

Qix = gGixi+g^Gix2 + -, (21) 

one obtains equations (for convenience of notation, the prime is now put 
outside the curly braces, which is justified because eigenvalues of Go are 
independent of A) 

u'xi = {(l-/)am4 ' (22) 

<2 = UXI {(1 - DGlXlYg, + {(1 - f)GlX2yg, , (23) 

etc., with solutions of the form [36] 

uxi = Kl-mxi}g, , (24) 
1 1 f"^ 

UX2 = 2"Ai + 2 y [usi,u',i] + {(1 - f)GiX2}go , (25) 

etc. This expansion can be rewritten in terms of the effective coupling gx for 
the purpose of eliminating ultraviolet divergences involved in the definition 
of the bare coupling constant and the corresponding counterterm. However, 
for g = gx~\~ Q(ff^), one can simply replace g by gx in the terms explicitly 



listed in Eqs. (24) and (25). In addition, Ux is unitary by construction. 
Therefore, Ux = Ux{boo) = Ux{bx) and one can freely replace 6oo by 6a in 
these equations, obtaining thus first two terms in expansion of Ux in powers 
of fi-A- 



^'^ This way of proceeding is not necessary, but it considerably simplifies the discussion 
that follows because there is no need to expand the form factor f\ and the energy 
denominators in powers of the coupling constant. 



14 
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4-2. Perturbative expansion for W 
Now consider Ac and A > Ac- Eq. ^ implies W{boo) = Ux^{boo)uj^( 



while in Eqs. (12) and (13) the operator W acts on states created from 
the vacuum by operators b\. It is therefore convenient to use an equivalent 
expression 

Wib^) = uiibx)UxM) (26) 

= 1 +5A'"ii +5A'"i2 + ••• ['^ + 9xu\,i+ glux,2 + ■■■] , (27) 

where 6oo in every u on the right-hand side is replaced by bx. Thus, 

W = l + gxWi+glW2 + (28) 
= {{fx-fxJGiooi}g, , (29) 



W2 = ^{uxi - ux.if + ^[ux,i,uxi] + ^ I ds[usi,u,-^] 



Ac 



+ {(1 - h.)QlX.2}g, - {(1 - fx)QlX2}g, , (30) 

where 600 is replaced by 6a everywhere on the right-hand sides. 

The perturbative expansion can be similarly carried out to higher orders. 
The term of focus here is Wi, as an illustration of how W acts on quark 
and gluon states. This illustration is valid when the coupling constant gx is 
very small, which means that A ^ ^qcd in the RGPEP scheme. 

To some extent, Wi is also useful as an indicator of how W acts on 
states when A approaches Ac ~ Aqcd- In this case, gx is relatively large. 
Although the size of contribution of terms containing Wk with A; > 1 in 
comparison to the size of contribution of the term with Wi for such small A 
is not known, it is certain that all these terms together contribute when 
A = Ac. They vanish no matter how large is gx because integrals in them 
effectively range from Ac to A and thus vanish when A tends to Ac- On the 
other hand, the higher power of gx the larger number of integrals involved. 
All these integrals must effectively range from Ac to A, so that they vanish 
when AA = A — Ac tends to as the appropriate power of AA. The first 
term in the expansion in AA must be linear in AA and Wi is such. 

Regarding the size of gx, the range of A right above Ac is a region where gx 
in the RGPEP scheme in QCD may be limited in size, instead of having large 
values that one may expect on the basis of a straightforward extrapolation of 
low-order expressions obtained in the region of large A [37] . Such limitation 
of gx from above is certainly observed in simple mathematical models that 
include asymptotic freedom and bound states in a renormalization group 
procedure for Hamiltonians [38j . Similar limitation of gx is also suggested in 
phenomenology of non-perturbative running couplings based on AdS/QCD 
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ideas and LF holography [39]. New theoretical information about the size 
of g\ at small A is expected to follow from a non-perturbative RGPEP 
equations described in Appendix [C| 



4-3. Wi for quarks and gluons 



According to Eq. (29), the first-order term in W 
is given by 



1 - {(/a - /aJ^/ooi}^(, , 



1 + gxWi + Oigl), 



(31) 



where TIq is the term independent of the bare coupling constant g and Hiooi 
is a regulated interaction term proportional to the first power of g in the 
canonical LF QCD Hamiltonian, Hoo = ^Ooo +5^/001 +0((7^), derived using 
gauge = from the Lagrangian for QCD, C 
Thus, 



^l){iIP-m)1p-\TrF^"'F|,y. 



Uoo = j ^^^^ [^Ooo + ghiooi + 0(3^)] , 
where the Hamiltonian density terms independent of g are 



(32) 



hoo 



1 



2 " ' id^ 

and the interaction terms order g are 



- TrA^{d^fA^ 



hi 



ool 



(33) 



(34) 



1 rab 
2" • 



with 5f/(3) color notation A = AH"", = if^H"", Tr{tH'') 

The calculation of Wi starts with the above expressions and proceeds as 
described in Appendix [B| 

Appendix [B] also argues that V7 for A > Ac ~ ^qcd contains all the 
terms in T-L\ that change the number of particles. This means that W 
can create a coherent slew of the effective Fock components at scale A in 
Eqs. (12) and (13) that are implied by the 2 or 3 constituent quark wave 
functions at Ac in Eqs. (10) and (11). In turn, this means that the result of 



action of W is fully encoded through the Hamiltonian in the wave functions 
■0Ac(l) 2) and ip\^{l, 2, 3) that correspond to the CQM. These wave functions 
are solutions to the eigenvalue problem for Hx^, which does not change the 
number of constituents. But the larger A the closer Hx to the canonical LF 
Hamiltonian of QCD and the more complex the states generated by W in 
Eqs. pl| and (pi). 



16 



SDGAPP2011GluonCondensate printed on July 1, 2011 



5. Eigenvalue problem for light hadrons 



According to previous sections, the same eigenvalue problem for light 
hadrons can be written in different ways when one uses different Hamil- 
tonian expressions that correspond to different values of A. Three values 
are distinguished in further discussion: A = oo, A > Ac, and A = Ac- The 
corresponding eigenvalue problems are equivalent if RGPEP equations are 
solved exactly. However, since exact solutions are not available, one is forced 
to guess plausible candidates for the Hamiltonian expressions at different 
values of A. 

At large A, insight comes from the perturbative expansion of RGPEP, 
using regulated canonical LF QCD Hamiltonian with counterterms to start 
with at A = oo and taking advantage of asymptotic freedom. At small A, 
one has to guess a first approximation. Phenomenology suggests that some 
form of a Hamiltonian for constituent quarks bound in a potential well is 
a good candidate for a first approximation to the QCD Hamiltonian with 
A = Ac for hadron states of smallest masses. 

These two extreme regions of A need to be connected to each other 
within the RGPEP framework in terms of some interpolation!^ The size of 
calculable corrections to any first approximation constructed this way will 
eventually tell us how far from a true RGPEP solution such first approxi- 
mation can be. The candidate for the first approximation that is developed 
in this article suggests that the concept of a gluon condensate may be in- 
terpreted in a new way. 



The first of the three ways, which uses the regulated canonical Hamil- 
tonian for LF QCD with all due counterterms, Hoq, amounts to formal 
writing of the eigenvalue problem in terms of Fock states created by prod- 
ucts of operators bto, dlo, and aL from the vacuum state |0) in the limit of 
regularization being removed, 



In this equation, a reasonably accurate description of the eigenstate 
presumably requires a very large number of wave functions for many sig- 
nificant bare-particle Fock components. The number of such components 
is expected to grow when the regularization is lifted. Although the full 
structure of is not known precisely, one can assume some model for one 

A similar reasoning applies also in the case of Yukawa theory, which is not asymptoti- 
cally free, when one considers a limited range of scales and effective coupling constant 
is sufficiently small from theoretical point of view to use perturbative RGPEP [40] . 



5.1. Three ways of writing the eigenvalue problem 




(35) 
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part of it, such as a Fock sector with a smallest possible number of bare 
particles in a meson or a baryon, and try to determine another part, such 
as a component with one additional gluon, using perturbation theory. Such 
approach is useful in description of exclusive processes that involve large 
momentum transfers [lljP^ 

The second of the three ways is designed for the opposite end of the 
scale for A in RGPEP, i.e., when A = Ac. The same Hamiltonian is expected 
to be expressible there in terms of operators for constituent quarks. In this 
case, one has the eigenvalue equation of the form 

HxM') = Ph\i^)^ (36) 



in which the state of a single light hadron is represented by Eq. (10) 



for mesons and Eq. (11) for baryons. This means that the light hadron 
eigenstates are described by only one wave function, ^/'c(l...n), with n = 2 
or n = 3. These wave functions appear in the Fock states built using 
creation operators h\ and . There are no components with effective 
gluons or additional quark-anti-quark pairs. In other words, the complexity 
of structure of light hadrons in QCD is encapsulated in the structure of 



constituent quarks. This means that the operators and in Eqs. (|10|) 



and ( |ll[ ) are complex combinations of products of operators 6^0, doz 
and their conjugates. 

Using Eqs. ([s]), (10) and (11), the eigenvalue problem of Eq. (36) for 
the n = 2 or n = 3 constituents is obtained in the form 

{^^Pi + ycn/pi^ xAl-.n\^) = P,-A.(l...n|^). (37) 

The unknown element in this form is the interaction potential term Yen- 
In order to derive the structure of Vcn in mesons and baryons, we shall 
employ the third way of writing the same eigenvalue problem using RGPEP. 
The third way suggests the reinterpretation of gluon condensate that is 
the subject of this article. When all gluons considered in the third way 
are incorporated in the constituent quarks at Ac, one is left only with the 
potential Vcn, see Appendix [D| 

A different way from RGPEP to attack the eigenvalue problem of a LF Hamiltonian 
for QCD head on is to use discretized light-cone quantization (DLCQ) in which one 
introduces a finite minimal unit of momentum, e^, which is the inverse of the size 
of a periodicity box for fields as functions of the position variable x~ on the LF. Since 
the total P"*" of an eigenstate of Hoc is conserved, the unit e+ naturally limits the 
number of bare particles in an eigenstate from above by the ratio P^/e^ [421 1431 HI] . 
Additional cutoffs need to be imposed in DLCQ in order to limit transverse momenta 
of constituents that may appear in an eigenvalue problem with some simultaneously 
fixed eigenvalues of the total momentum components and P^ . 
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The third way of writing the same eigenvalue problem uses Hx with 
A > Ac- The size of A is assumed to be such that gluons are already some- 
what active as participants in the dynamics while the effective quarks are 
still heavy enough for a NR analysis to apply to their relative motion in 
a light hadron. The NR approximation is thinkable for a slowly moving 
hadron according to previous Sections, since interactions in the effective 
Hamiltonian at scale A > Ac contain form factors. The effective particles 
of sizable masses cannot change their momenta through interactions by 
large amounts. Thus, an eigenvector corresponding to a slowly moving light 
hadron can be expected to be described by the wave functions that have 
significant values only for small momenta of the constituents. 

The third picture at some scale A > Ac is considered a candidate for a 
first approximation to hadrons in RGPEP. The corresponding Hamiltonian 
will be proposed below using ideas motivated by gauge symmetry in its NR 
form. Difference between this approximate NR form and the interactions 
that can be systematically calculated using RGPEP in LF QCD, is to be 
treated as a perturbation^^ The approximate picture at A > Ac is described 
in more detail in the next Sections. Here we only suggest that RGPEP 
provides a scheme that may be used in future calculations of Vcn in Eq. 
(37) at Ac by evolving the third- way picture at A > Ac down to Ac and 
taking advantage of the RGPEP reinterpretation of the gluon condensate 
that is described in the next Sections. 



5.2. Light hadron states at X> Xc 



When A > Ac, the meson and baryon states are described by Eqs. (12) 
and (13), respectively. In these equations, the operator W acts on the Fock 
components with 2 or 3 quarks of scale A and creates additional components. 
The meson state is simpler than the baryon state and it will be discussed 
first. 



This strategy is similar to the one described in Ref. [iGi with an artificial linear po- 
tential. RGPEP introduces new elements: the dynamical transformation from bare to 
effective particles available at different scales, construction of corresponding effective 
fields, possibility of using NR approximation at A > Ac for identifying interactions 
through gauge symmetry as indicated later in the text (instead of introducing an 
artificial potential), extension beyond perturbation theory (see Appendix [c|, and 
invariance with respect to 7 kinematical LF symmetries for the resulting interaction 
terms and wave functions [27]. As a result, the first approximation potential is pro- 
posed to be not a linear but a quadratic function of a relative distance between color 
charges (see next Sections). The quadratic potential is expected to cause creation of 
additional particles when a distance between two colored particles increases and the 
net linear increase of energy with distance must ultimately result from the energy of 
particles created along a line that connects the two widely separated color charges. 
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The meson state is 



(38) 



12 



From Eqs. (B.16), (B.18) and (B.22), it is clear that W shares properties 
of the interaction Hamiltonian. For example, if Hxj changes a group of 
particles to another group, a similar change with additional factors results 
from action of W . So, out of 2 effective quarks in a meson (3 in a nucleon) 
more quarks and gluons are created. In order to use local gauge symmetry 
to propose the approximate form of the state that results from action of W , 



it is useful to represent the state of Eq. (38) in position space. This is done 



using quantum fields built from effective particle operators at scale A. 



The effective quantum fields are constructed using Eq. (B.l) for quarks 
and Eq. (B.8) for gluons and replacing operators h, d and a with h\, d\ and 



ax, respectively. Thus, the operator ipxix) on the LF is built in the same 



way from bx and d^^ as the canonical operator ■il}{x) is built in Eq. (B.l) 
from the operators b and d^ that are equal boo and dto, respectively; 



ipx{x) 



E 

acf 



[k] XcUXfkabkacfXe ^''^ + XcVXfkud\^^jy^ 



Akx 



(39) 



The spinors uxfka and vxfka include corresponding vectors in flavor space. 
Spinors might depend on A if the effective quark masses in the field expan- 
sion are allowed to depend on Aj^ Similarly, operator A^{x) on the LF is 

built from ax and as the canonical operator A^(x) is built in Eq. (B.8) 

from the operators a and a^' that are equal a^o and aL , respectively. Gluon 
polarization vectors do not depend on A. 



ac 



[k] 



4-C ^ 



-ikx _|_ J.C t ikx 



(40) 



As a consequence, the dynamically independent components, ipx = 

and A-^ in = gauge on the LF = have the same commutation 

relations as in a canonical theory. 

With particle operators at A = Ac, an effective constituent quark field 
operator is constructed in the same way. When this field is used to describe 
a slowly moving meson, it is useful to write the field as 



Ac 



(41) 



^* Note that the independent field components, = A^-)/), A"*" = (l/2)7''7"'", are 
independent of the quark masses, and inclusion of effective masses in spinors is merely 
a way of useful notation for some effects of interactions. 
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where the upper two-component field U\r annihilates quarks and the lower 
two-component creates anti-quarksrj 

Thus, by inverting the Fourier transforms through integration over the 
LF hyperplane and using conventions explained in Appendix [B| the state 



in Eq. ( 10 ) for a meson of momentum Ph can be written as 

X e-'(*^i^i+'=^^^)C/l (xi)V2x2(M^A.(x2)|0), (42) 



where V'2x2(fei2) denotes the 2x2 matrix wave function of relative motion 
of the quarks. The three-vector ki2 can be defined as a relative momentum 
of the quarks in the CRF. Details of the definition of ki2 are not essential 
at this point but later discussion will include relevant details. 

In order to obtain analogous position representation of the meson state 



in Eqs. (12) or (38), one needs to apply W to the expression in Eq. (42). 
The result is 

\^)m = ^ / ^^^^ / [12] lQn'P^5\P, -k,- k,) 

X e-*('=^^^+'=^^^)VFC/l(xi) ^2x2(^12) 1^a(^2)|0). (43) 

A similar expression is generated in terms of three quark fields and W for 
baryons. Integration over momenta in these expressions generates position 
wave functions as coefficients in the expansion of meson or baryon states 
into basis states that are created from the LF vacuum by action of a product 
of two or three quark fields at scale A and W . 

Continuing with the meson states, the central hypothesis about action 
of W is that gauge symmetry forces the result of action of on a state of 
two quarks at scale A to have the form 

W J d^xid^X2Ul{xi)i:{xi,X2)Vx{x2m (44) 

d^xid'^X2d^xo Ul{xi) W{xi,X2, xq) Vx{x2)\0) , 

(45) 



The three- vector notation x or k refers here to the LF co-ordinates (a;~ , a;^ ) in position 
space. A similar notation is sometimes also used for momentum variables k^). 
However, when we proceed later to the Schrodinger eigenvalue problem for effective 
Hamiltonians for light hadrons, the same notation will be adopted also for three- 
vectors built from _L and -I- or — components in such a way that the standard three- 
dimensional notation respecting rotational symmetry will be natural. 
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where the operator W{xi,X2,xq) is responsible for creating components 
generated by W . These components can be of two kinds. One kind is 
formed by operators that carry the color of initial two quarks. For example, 
if W generates a gluon from a quark, the generated gluon and the emerging 
quark together carry the color of the initial quark. The other kind is formed 
by colorless operators. For example, W may cause emission of two gluons 
by a quark and the gluons may form a color singlet. 

These two kinds of contributions are encapsulated in W{xi,X2^xo) by 
writing 

W{xi,X2,Xq) = ^{xi,X2,Xo)T{xi,X2)G\xq) , (46) 

where iI:{xi,X2,xq) denotes a new wave function at scale A > Ac, G'^{xq) 
denotes the colorless component of the state, and 



r(xi,X2) = -Pexp 



-ig / dx^A^(x) 

J X2 



(47) 



is the color-transport factor along a straight line between quarks that main- 



tains local gauge symmetry by bringing in required gluon fields ^® This 
factor is constructed in analogy with Ref. [17j. 

The operator that generates the colorless component of a hadron state, 
(jr^(xo), will provide the contribution in dynamics of quarks that is associ- 
ated with gluon condensation in hadrons, rather than in a vacuum. Namely, 
instead of the vacuum expectation values of operators considered in Ref. 
[T7] . such as (ril^d^A-^ with i,j = 1,2, the Schrodinger equation for the 
wave function ip{xi,X2,xo) will involve expectation values 

|G> = Gt(xo)lO). (49) 

A similar reasoning is followed regarding baryons. In the case of baryons, 
one has to deal with three color-transport factors that are constructed in 



In the LF gauge A'^ = 0, one may expect only transverse separation between quarks 
to count. However, the dependent (constrained) components of fields, ip- and A~ , 
contribute to the effective interactions in a non-trivial way and one has to keep in mind 
that the complete effective theory at small A should have full rotational symmetry 
restored. Therefore, a complete RGPEP expression for W must account for the 
dynamics along — as well as along xj^ — X2 ■ This issue will be addressed later 
by introducing effective interactions that respect rotational symmetry inside slowly 
moving hadrons through a new definition of three-dimensional relative momenta of 
constituents to which the minimal gauge coupling rule can be applied in a rotationally 
symmetric way. 
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analogy to Ref. |17j . The operator in baryons generates the state G^|0) 
that plays the same role that the vacuum state jO) played in Ref. jl7j . 
It is assumed that the colorless components of mesons and baryons are 
approximately the same. The issue of universality of expectation values 



such as in Eq. (48) for mesons will be further discussed below when we 
come to the construction of the effective Hamiltonian. 

In summary, the claim of gauge symmetry regarding the colorless basis 
states of quarks and gluons at scale A from which mesons and baryons are 
made, is that they are of the form (all fields are effective at scale A > Ac) 

s:ab 

\xuX2,xg) = ^^{u\^nr{T^V2xfG^\0), (50) 



ab ^ 



^abc 



\Xi,X2,X3,Xg) 



where 



^ (ulT^r {u\T,f {ulnr Gt|0) , (51) 

abc 



^-igQdx.A^^ (52) 



and X = (xi + X2)/2 in a meson and x = (xi + X2 + X3)/3 in a baryon. 
The factor T is defined here along a straight path in such a way that it 



reproduces the color-transport factor in mesons in Eq. (47). The path 
dependence is ignored as an unnecessary complication at the level of mean- 
field approximation. The operator G with argument xq represents the white 
energy density background called glue. It is assumed to be a scalar boson 
field with corresponding commutation relations. 

There is a trouble on the LF with non-locality of the boson commuta- 
tion relations in the direction of x~ , which requires a solution. However, 
the effective dynamics that will be constructed in next Sections will explic- 
itly circumvent this difficulty by using an operator G that is a function of 
a three-dimensional position space variables associated with slowly moving 
hadronic constituents in a slowly moving hadron. Thus, the problem ulti- 
mately requiring a solution is not so much the non-locality of a scalar field 
but how LF QCD can generate a rotationally symmetric effective theory for 
light hadrons. The construction offered in next Sections proposes to treat 
G as a field depending on a suitable three-dimensional position variables in 
which it can be local in a way that respects rotational symmetry. 

The quanta of G are meant to represent excitations of the states of 
gluons condensed inside hadrons. The quanta of G are not point-like. Their 
size is characterized by 1/A and can be considered roughly on the order of 
1/Aqcd, i-e., the quantum extends over the volume of an entire hadron 

^"^ The effective glue degree of freedom represents contributions of all white Fock sectors 
of effective particles that share the hadron momentum to a varying degree as A varies. 
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There exists a possibility to associate vector or even higher spin nature with 
G, contributing to the hadron spin, but there is no need to do so here. 

The basis states are orthonormal in the sense that, for all quarks being 
different, 

{xi,X2,Xg\xi',X2i,Xg') = '5ll'(^22"^GG' , (53) 

{xi,X2,xz,XG\xv,x2i,xy,XG') = hv ^22' ^zy ^GG' , (54) 



and 5kk' includes S'^{xk — affc')]^^ If quark quantum numbers besides color 



are not different (this comment concerns only the baryon case), one has to 
adjust the normalization of basis states by including all permutations that 
contribute to the scalar products. 

In the abbreviated notation used in the remaining part of the article, 

\xi,X2,xg) = |12G), (55) 
|^i,^2,^3,^g) = |123G). (56) 

In these abbreviated notation, the meson and baryon states read: 

Y,^{UG)\12G) , (57) 

12G 

J]] ?/;(123G) |123G) . (58) 



IM 

12G 



)B = 

123G 



5.3. Preliminaries concerning Hamiltonian density at X> Xc 

Using the concept of effective quark and gluon fields, one can propose 
that the LF Hamiltonian at A > Ac be expressible, by analogy to the canon- 
ical theory, in terms of an integral over the LF of a density that is a function 
of the fields. Thus, 

Hx = f {hox +hix) , (59) 
Jlf 

where the density /iqa is bilinear in the fields and renders single-particle 
operators while the density hjx describes interactions in terms of products of 
at least three fields (the interactions involve at least three effective particles) . 

Using power-counting [16] , the bilinear density at scale A can be assumed 
in the form (all fields at x"*" = 0) 

hox = li^xl^ ""V a + TrAii-d^ ' + m^)^f + CTix , (60) 

With the quahfication that locahty in the x~ direction requires proper definition 
of z-components in a rotationally invariant efTective theory to be discussed in next 
Sections. 
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where the masses correspond to A. The symbol CTix denotes mass-like 
terms needed to make the mass parameters m and nig to count as masses in 
the eigenvalue equations for light mesons and baryons. In other words, CTjx 
are by construction cancelled by self-interactions in light colorless states 
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The interaction density hjx in Eq. (59) must be non-local in the sense 
that it involves products of fields at different points In particular, the 
minimal coupling between quarks and gluons, which in the canonical gauge 
theory is of the form 



J d^xhMCoo{x) , (61) 



HmCoo 

with 

hucooix) = gi){x) 4{x)%l){x) , (62) 

in the effective theory must take a non-local form that in a lowest-order 
approximation is of the type |27| 

Hmcx = j d^xi Sx2 d^X3 hMC\{xi,X2, xs) , (63) 

with 

hMC\{xi,X2,X3) = g\f\{x2 - Xi,X3 - Xi)i)\{xi)4.\{x2)'ip\{x3) . (64) 

This non-local interaction term appears with other non-local terms in the 
Hamiltonian obtained from RGPEP at scale A. The dependent and inde- 



pendent components of the fields are grouped in Eq. (64) according to a 
free theory. The dependent components involve the inverse of id'^ , which is 
a non-local operator. This non-locality is of the same type as in a canonical 
theory and the RGPEP non-locality appears here on top of the canonical 
one. 

The non-local interaction terms with small A can be simplified consid- 
erably if the domain of action of the Hamiltonian is restricted to slowly 



There is no need to specify these terms further here. Such terms can be identified 
exphcitly in perturbation theory, e.g., proceeding in a similar way to how it is done 
in the eigenvalue problem for heavy quarkonia in Ref. |18| . 

The non-locality discussed here is due to the RGPEP vertex form factors of width A 
in momentum space. The RGPEP-induced non-locality must not be confused with 
the canonical non-locality of LF Hamiltonians in which the dependent parts of fields, 
such as ij)- = K-tj) and A~ , involve inverse powers of id'^ , which is a non-local 
integral operator. 
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moving effective quarks. Namely, consider again the case of Eqs. (63) and 



( 64 ) and introduce a gradient expansion of the form 



hMCx{xi,X2,Xz) = g\f xix2 - Xi,X3 - Xi)i!x{xi) 

X [Mxi) + dMxi){x2-xi) + ...] 

X [^xixi) + d'ipx{xi)ix3 - xi) + ...] . (65) 

The three dots indicate terms with higher derivatives. If all effective par- 
ticles move slowly and the derivative terms are small, the Hamiltonian can 
be approximated by the first term in the expansion, 



Hix = J d'xi d^X2 d''x3 hix (66) 
^ gx [ d'^xMx)Mx)Mx) I d'^yd'^zfxiy,z) + .... (67) 



This means that the non-local effective interaction in a slowly moving, NR 
system still looks like a local one except that its strength is determined not 
solely by the coupling constant gx but also by the integral of the non-local 
form factor on the LF hyperplane, 

gx = gxj d^yd^zfx{y,z). (68) 

The point is that the effective Hamiltonian density at A > Ac may 
partly resemble a Hamiltonian of local gauge theory in its terms that couple 
quarks with gluons even though the actual effective interactions are non- 
local, provided that one limits the domain of the Hamiltonian to slowly 
moving hadrons. In this case, construction of approximate candidates for 
the LF Hamiltonian density of coupling between quarks and gluons, before 
they are corrected using RGPEP, may proceed in analogy to QED |47] . This 
means that one uses the gauge j4+ = 0. The derivative id-^ in the quark 



kinetic energy in Eq. (60) is supplied with the minimal coupling addition 
of gA-^ (the product of derivatives is separated by the inverse of id~^ , which 
is not altered). In addition, one includes terms dictated by constraints that 



imply the result of Eq. (62) for the quark-gluon coupling. 

Much less is understood about the gluon part of the effective Hamil- 
tonian density for light hadrons at A > Ac. The lack of understanding of 
the Hamiltonian is reflected in the lack of understanding of the gluon com- 
ponents in its eigenstates. One would have to calculate the operator Hx 
precisely in order to uncover the information it contains about how to build 
a model approximating light hadrons in LF QCD. While RGPEP provides 
tools for such calculations, the remaining part of the paper is only devoted 
to deriving the model that can be treated as a first approximation. 
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A simple candidate for an approximate Hamiltonian for light hadrons 
in LF QCD is constructed in the next Section assuming that: (1) the result 
of action of W in RGPEP can be represented by inclusion of the color- 
transport factors T that maintain local gauge symmetry of hadronic states, 

(2) the operator G represents the condensation of gluons inside hadrons, 

(3) a mean field approximation can be applied to the gluon field operator 
A in the minimal coupling of effective quarks to the gluons condensed in a 
hadron, and (4) the effects of quark and gluon binding that are not treatable 
in the mean field approximation can be included by an ad hoc Gaussian 
approximation to the wave function of relative motion of constituent quarks 
with respect to the condensed gluons. 



5.4- Mass squared with minimal coupling at X> Xc 

The leading idea is that the approximate LF Hamiltonian for 2 or 3 
quarks at scale A > Ac in light hadrons should have a form compatible with 
minimal coupling between quarks and gluons and it should respect Poincare 
symmetry. Realization of the idea employs the four assumptions listed at 
the end of the previous Section in the following way. 

The LF Hamiltonian eigenvalue problem at A > Ac for a light hadron 
built from 2 or 3 quarks and condensed gluons, 

HxH) = ^^IV'), (69) 

can be written in terms of the eigenvalue equation for the effective invariant 
mass operator of the hadron, M\, in the form 

Mll^l^) = M^\^). (70) 



One can think about the operator on the left-hand side of this equation in 
terms of the invariant mass squared of free particles plus interaction terms. 
For a free quark of mass mi and a free anti-quark of mass m2 , one has 

M%ee = {\]ml+pl + sjml+pi^ - {pi +P2f . (71) 

For slowly moving particles, neglecting terms smaller than the ones that are 
quadratic in momenta, the NR approximation renders 

Mfree = {mi + rUi) + — " Pi + " P2 - Pqq , (72) 
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where P denotes the total momentum of the two quarks. Using 

/3i = ^^, (73) 



mi + m2 

h = (74) 

m\ + 7712 



one has 



-M/ree = ("^1 + "^2)' + ^ + ^ - . (75) 



According to the gauge rule of minimal coupling, the interaction of quarks 
with the condensed gluons should be described by 

V\ - 9^\\ (p2 - 9^2^ 



^minimal = ("^1 + ^2)'^ + ^ 5 — + ^ o ' ^q] > (76) 

Pi P2 

where Ai is an abbreviation for A{xi). Analogous reasoning in the case of 
baryons yields 

/ 3 \ 2 3 f^,_g^,Y 

where /3j = mj/(mi + m2 + rn^). 

For all quarks having the same mass m, (3m = /3i = /32 = 1/2 in mesons 
and /3_B = /?! = /32 = /^a = 1/3 in baryons. Assuming these simplifications, 
the minimal coupling terms in Eqs. (76) and (77) can be interpreted follow- 
ing Ref. [17j in the context of standard Hamiltonian dynamics as resulting 
from the Hamiltonian operator with a proper (3, i.e., /3m in mesons and /3b 
in baryons, 

Hmin = d^x : V;(x)HV^. -5A(f)]V(^) : . (78) 

5.5. Reinterpretation of the gluon condensate 
Using Eqs. (|57|), (1581), and (ITSI), one can consider the eigenvalue problems 



H^inY,^{12G) \12G) = ^^^(12^) \12G) , (79) 

12G 12G 

Hmin ^(123G) I123G) = M|, ^ %b{l23G) |123G) , (80) 



123G 123G 
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for quark subsystems in mesons and baryons and project them on the cor- 
responding basis states. The calculation proceeds in a similar way to the 
one in Ref. [17], with three exceptions. 

The first difference is that the eigenvalues M^^ and Mg^ refer only to the 
qq subsystem in mesons, instead of the entire meson mass squared, and to 
the 3q subsystem in a baryon, instead of the entire baryon mass squared. 
Note also that the eigenvalue equations considered in Ref. [T7] were for 
usual energies while we now consider the invariant mass squared operators 
that include the quark free invariant mass squared and minimal coupling 
interactions. The second difference is that the vacuum state \^}) is replaced 
by the state of gluons condensed in a hadron, 

\G) = Gt(fG)lO). (81) 



The vacuum state |0) does not develop any condensate. The third difference 
is that in the mean-field approximation for the gluon field 



21 



A{x) = ^B{xg) X (x-xg), (82) 

one introduces the color magnetic field operator at the center of the gluon 
component G of a hadron, B{xg), instead of the vacuum operator B at an 
arbitrary point xq. However, these differences do not change the formal cal- 
culation of the matrix elements. The field B{xg) is assumed factorized into 
a product of vectors, one with space and the other with color components 
as in Ref. [IT] . This assumption is meant to reflect the assumed lack of cor- 
relation between position space and color space directions in the mean-field 
approximation and it renders a simple, Abelian form of the approximate 
model for the effective Hamiltonian. Matrix elements of terms linear in A 
are set to zero because they change under gauge transformations whereas 
the white component of a hadron does not. 

The matrix elements one obtains for a qq pair in a meson and 2>q in a 
baryon reacj^ 



{l2G\H^in\ll>)M 





This approximation replaces the Schwinger gauge formula for the background gluon 
field that reproduces QCD sum rules for heavy quarkonia in the LF Hamiltonian 
formulation of the theory |48| . 
■^■^ Eq. (14) in Ref. [T7] misses the factor 1/3 in front of the condensate term that is 



correctly printed here in Eq. ( 83 \ 
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1 1 3 2 

{12'iG\H^MB = V(-A,)(123G|V^)b + — 

X Tr^ ^-^^ (123G|^).. 

(84) 

All gluon field matrix elements contain squares of differences of the effective 
gluon field operator at different points. Therefore, the position of the gluon 
body, xg, drops out. In the mean-field approximation, one has 

A{x)-A{y) = is X (x-y), (85) 
where B = B{xg)- Consequently, 

Tr{G\^{A,- Ayf\G') = 2{x - yf Cgiue{G\G') , (86) 

with 

{G\G') = 6gg' . (87) 

The expectation value Cgiue plays here the same role that the vacuum gluon 
condensate value Cyacuum = ^lacuuml^^ with 

^lacuun. = (f^|(a/vr)G^-G^jn), (88) 

plays in Ref. [17] . Replacement of the constant Cyacuum by the constant 
Cgiue = ^1iue/^^-> implies our reinterpretation of the phenomenologically 
useful quantity vacuum on the order of Aqcd as coming from the quantity 
fgiue that originates according to the RGPEP model in the gluon conden- 
sation only inside a hadron instead of the entire space. Hence, 



{12G\Hmin\i')M = — (-Ai - A2)(12G|V')m 

2 



{123G\Hrmn\i^)B 



PM 


-Ai-. 


J-( 


'T^^glue 




. 3 






1 5 




^8 


V 3 



)'^(12G|^)m, (89) 



12 , ^'3 



(123G|^)b, (90) 
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where 

ri2 = XI-X2, (91) 
rs = x^-{xi + X2)/2. (92) 

The results for Hmin with new interpretation of the gluon condensate inside 
hadrons are used in the next Section to construct candidates for effective 
LF Hamiltonians of hght mesons and baryons as first approximations to 
solutions of RGPEP in QCD. 



5. 6. Approximate LF Hamiltonian for quarks at X> Xc in light hadrons 
Candidates for approximate Hamiltonians for light hadrons can be ob- 



tained starting from inclusion of results in Eqs. (83) and (84), or (89) and 



(90), in Eqs. (76) and (77), correspondingly. The meson case is simpler than 



the baryon case and explains the leading idea of constructing the effective 
Hamiltonians using the minimal coupling dictated by gauge symmetry. We 
have 



M 



Using the variables 



+ 



^ \2 



{mi +1112) -{pi+m) 



+ 



1 
1 

J2 



3 



pI + ( '-Tr{Ai - A2'^ 



Vl + ^Tr(Ai - A^f 



P12 
k 

Pi 

P2 



PI+P2, 

/32P1 - P1P2 . 
/3iPi2 + k, 
hPi2-k, 



G. 



G. 



(93) 



(94) 
(95) 
(96) 
(97) 



one arrives at 



'■qq 



(mi + 7712)^ + 



/3i/32 



-Tr{Ai 



A2? 



G. 



(98) 



The LF counterpart of this result is obtained by considering three-vectors 
with -|- and _L components instead of z and _L. So, one writes 



12 



xi 



Pi 

Pt/P^ 



+ P2 , 
12 ' 2:2 



Pt/P^ 



12 ' 



(99) 
(100) 
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Pi 

vi 



X2P1 -X1P2 



12 



(101) 
(102) 
(103) 



These are standard expressions in LF description of relative motion of two 
constituents carrying the total momentum P. Using these variables for two 
free particles named 1 and 2, one obtains their free invariant mass squared 
in the form 



12 



+ — 



{mi + m2)^ + 



X2 



X1X2 



+ (m2Xi - 7711X2) 



(104) 
(105) 



Comparison with Eq. (98) identifies the relative momentum k, 



k' 



X1X2 
XIX2 



[7772X1 — 7771X2] 



(106) 
(107) 



This is a new way of parameterizing relative motion of two constituents 
in LF approach to quantum mechanics and field theory. The new vari- 
ables match non-relativistic relative momenta that include effects of parti- 
cle masses. On the other hand, it is known that in the relativistic relative 
motion of two constituents the mass parameters appear not significant. For 
example, the standard relative momentum, in which fc"*" = k'^, has length 



k^ 

^standard 



4Mi2 



[MI2 



(7771 + "72) 



[MI2 



(7771 - "72)2] , (108) 



and a relatively complicated expression is obtained for k^ if one insists on the 
identification k-^ = k-^. But when A^i2 is large, one has k'^tandard ~ -^12/^ 
independently of the quark masses. At the same time, the variable k^^^^^^^^ 
depends on k"*". The length of the standard relative three-momentum in the 
constituent rest frame differs from the length of the new one. 



standard 



4/3i/32 



(7771 — 7772) 



21 



(109) 



and the angular orientations of the three-vectors kgtandard and k also dif- 
fer. This means that the new choice of LF three- momentum variables in 
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relativistic cases involves rotation by some polar angle (the azimuthal an- 
gles are the same). Analysis of rotational symmetry is influenced in the 
sense that in order to recover a simple NR quantum mechanical picture of 
hadrons from quantum field theory one is motivated by the effective picture 
in RGPEP to use the new variable k rather than kstandard as arguments 
of potentials for constituent quarks. One should also remember that when 
a change of variables is made in the expressions involving effective quark 
wave functions, the new variables produce Jacobian factors in phase-space 
integration. 



Using Eq. (89), one obtains 



9 



Tr{Ai - A^)' 



G 



1 ( T^^glue 

2 V 3 



12 



(110) 



where r\2 should be the relative position variable that is canonically conju- 
gated to k. This means in quantum mechanics that 



r\2 



. d 
dk 



(111) 



Thus, the LF mass squared for a constituent quark-anti-quark pair at scale 
A > Ac interacting in a gauge minimal way with gluons condensed inside a 
meson and treated in a mean-field approximation, has the form 



M 



(mi + + 



1 



/3i/32 



1 / ""'(fglue 

2 V 3 



2 1 
2 



d 

dk 



(112) 



where the vector k is defined by Eqs. (106) and (107) 



In particular, the effective quarks u and d are expected to have prac- 
tically the same mass m order Aqqd in the RGPEP scheme. For them, 
/^i = = 1/2 and the associated Jacobi variables are 



P 

Pp 



n2/V2, 
kV2. 



In terms of these variables. 



■""■qq 



W + 2p^ + 2(!^, p- 



(113) 
(114) 



(115) 



For small relative momenta and corresponding distances, this result approx- 
imately matches the NR oscillator dynamics. 



M 



2m + 



2m 



+ 



1 



m 



V 3m 



(116) 
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with frequency ujm — T^fgiue 

/(3m). 

On the one hand, this result relates the RGPEP reasoning to the CQM 
phenomenology since the reinterpreted value of the gluon condensate pro- 
duces physically reasonable frequency ojm [E]- On the other hand, the 
variables k and r identified here include the characteristic factor \/x{l — x) 
that is required in AdS/QCD holographic variables [Ml |50l |51] • Thus, it be- 
comes plausible that also the RGPEP reasoning regarding the A-dependence 
of effective interactions may provide insight into the significance of soft wall 
(SW) models [52j of hadron spectrum. 

As a result of analogous reasoning in the baryon case, we obtain 

Ml^^-im^+m^ + m^f = ^ Q^ + 1^k\ (117) 

P3(,l-P3j P1P2 



where 



A= ^""V , (119) 
nil + 1TT.2 + m3 

pi = x^Phz + 9^ , (120) 

P2 = X2Ph3--^q^-k^, (121) 



123 + 


q~ 


L 






X2 


'123 ~ 


1 


- X'i 






Xl 


'123 ~ 


1 


- X'i 



Pt = ^1^23 - q^ + k^, (122) 



y 0:3(1 - X3) 



= J ^4^3^ [("^1 + "^2)^3 - m^ixi + X2)] , (124) 
y X^^L X3) 



K-=M^\-fk\ (125) 

V X1X2 I- P3 



//3i/?2l-X3 f Xl X2 \ , . 

For slow relative motion of constituents, the new momentum variables K 
and Q match the non-relativistic three-momenta used in quark models, but 
in fact they apply in the entire range of relativistic kinematics of motion of 
constituents inside baryons, and for arbitrary motion of baryons as a whole, 
thanks to the kinematical symmetries of LF formulation of the theory. How- 
ever, one needs to remember that the lengths and angular orientations of 
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the new momentum variables differ from the standard LF variables among 
which, in particular, and (see below) could be thought directly use- 
able for obtaining some effective constituent picture for baryons from LF 
QCD. 

The LF mass squared for 3 constituent quarks of one and the same mass 
m at scale A > Ac interacting in a gauge minimal way with gluons condensed 
inside a baryon in a mean-field approximation, has the form 

Ml = 9m' + + Iq' - 3m' liAl/2 + 2AQ/3) .(127) 

We identify relations 

K = qi2 = Pp/V2, (128) 

= ri2 = V2p, (129) 

oK 

Q = q3 = -VV^PX, (130) 

-i^ = rz = -V372A, (131) 

in terms of the Jacobi co-ordinates for 3 quarks with equal masses, p, pp, A, 
and p\. One has 

-Pl23 = P1+P2+ P3 , qi2 = {Pl - P2)/2 , (132) 
Pl2 = P1+P2, 93 = (2P3 - Pl2)/3 , (133) 

Pi2 = 2P123/3 - 93 , P3 = P123/3 + qs , (134) 

Pi = P12/2 + qi2 = -P123/3 + 912 - 93/2 , (135) 

P2 = -P12/2 - 912 = -P123/3 - 912 - 93/2 , (136) 

P3 = P/3 + 93. (137) 

pI+pI+pI = ^^123/3 + 29?2 + 39i/2 . (138) 

R123 = {xi + X2 + X3)/3 , ri2 = xi - X2 , r3 = X3- R12 , (139) 

X3 = i?i23 + 2r3/3, R12 = ixi+X2)/2 = i?i23-r3/3, (140) 

XI = R12 + ri2/2 = i?i23 - r3/3 + ri2/2 , (141) 

X2 = R12 - ri2/2 = -R123 - r3/3 - ri2/2 . (142) 



^PjXj = P123-R12 + 912^12 + 93?'3 • (143) 
i=l 
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In the case of the same masses, the relativistic LF result for three quarks in 
a bar yon reads, 



Ml^ = 9m2 + + + (|^)' ^(p^ + X^) . (144) 



In the NR hmit. 



which matches precisely the result of Ref . [T7] , and 

ul = ^ul,. (146) 

This numerical relation is in a reasonable agreement with phenomenology 
of constituent quark models [211 1221 [23] , as noticed already in Ref. [T7j (see 
also Section 6.4). 



5. 7. Dynamics of quarks and glue in hadrons at X> \c 

The constituent models do not, however, include the gluon component 
G. Model potentials are also not assigned energy or momentum, while glu- 
ons do carry energy and momentum. Bag models do include a bag in terms 
of boundary conditions on quark wave functions on the bag walls and a 
constant contribution of the bag to the total hadron energy, but the bag is 
not placed yet in the context of QCD [53]. The author does not know how 
to interpret results of lattice formulation of QCD regarding the constituent 
model and the gluon component in terms of hadronic wave functions. Sum- 
marizing, it is not clear how the addition of G will affect phenomenology 
familiar through the CQMs without G, although some foreseeable changes 
appear welcome (see next Section). 

Regarding the contribution of mass of G to the mass squared of the whole 
hadron, one cannot say anything but hypothesize that the mass of G may 
depend on A. It can be small, due to interactions that hold gluons among 
themselves in the form of G. The mass of G may even decrease when A drops 
down to Ac. A graphical picture that makes this idea plausible as a result of 
overlapping of quarks of size 1/A (a quark and an anti-quark in a meson or 
three quarks in a baryon) and G in position space, is described in Appendix 
[D| However, a simple concept of a constituent G may be insufficient in the 
case of vr-mesons, considering that vr- mesons may differ from other hadrons 
as a consequence of their relationship to chiral symmetry and its breaking. 
The symmetry breaking can be addressed in LF QCD [16] but it is not 
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addressed in this article. Such discussion requires reinterpretation of the 
quark condensate and G built from gluons alone is not sufficient. 

Regarding the motion of the glue component G inside a hadron, one may 
observe that if the state G^lO) in RGPEP is to correspond to the vacuum in 
the instant dynamics way of thinking and G is associated through RGPEP 
with a hadron in an effective theory of scale A > Ac, the effective quarks at 
this scale should not be free to move farther away from xq than about 1/A 
(see also Appendix [D]) . The mean field approximation misses the fact that 
the gluon component has a size comparable with a hadron. The Schwinger 
gauge expression for the gluon field potential Ax{x) with only one term 
linear m. x — xq and B{xg), Eq. ( |82[ ), must be missing important effects 
at distances comparable with the diameter of a hadron, cf. [48|. This 
means that the gauge transport factors T{x,x) cannot be approximated 
well in the entire volume of a hadron using the Schwinger gauge with a 
linear term alone. Thus, the mean-field calculation does not report well on 
what happens at the boundary of the glue or the boundary of a hadron. 
While one cannot determine in advance what will eventually result from 
RGPEP calculations of effective Hamiltonians and eigenvalue problems for 
these Hamiltonians, one should certainly expect that there will be forces 
that keep the center of the quarks and the glue together. 

There is no simpler choice for the needed interaction term between the 
quarks and G than a quadratic function of the distance between x and xq- 
Such distance is easy to write in NR quantum mechanics. However, one has 
to carefully define its analog in the hadron mass squared operator on the LF. 
Solution of this technical issue described here, is found by looking at hadrons 
as built only from two constituents. One constituent is the effective quarks 
treated as one particle of mass resulting from eigenvalues of -Mqg in mesons 
or M'^g in baryons. The other constituent is G treated as a particle. The 
mass of G is unknown, may depend on A, and is free to adjust in a process 
of defining a first approximation to hadrons at A > Ac. With this setup, 
one applies the same method for constructing the LF mass squared operator 
that we used for a quark and an anti-quark treated as two constituents in 
a harmonically bound state. There is no minimal gauge coupling to use for 
the colorless state of quarks interacting with a white G, and a large degree 
of ambiguity must be confronted without clear guidance from QCPj^But 
there is no problem now with introducing the appropriate oscillator potential 
operator using the distance that is quantum-mechanically conjugated to a 
relative momentum. 

The first-approximation to LF Hamiltonian for light hadrons at A > Ac 



It is precisely this type of ambiguity at small energies that RGPEP is meant to 
eventually help in resolving. 
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is suggested to have the form 

TT •'^ quarks^ -^quarks , •'^G'^^G i 9G /i /i'7\ 

^A>Ac = ^+ + ^+ + ^ ■ (147) 

quarks G hadron 

In order to define a suitable candidate for the last term, one can introduce 

^hadron ^ ^^uarks + "^g' ' i^^^) 

^1 ~ ^quarks/ -^hadron ' (149) 

^G = P^/PLdron^ (150) 

= XGPquarks ~ ^qPc ' (1^1) 

Pquarks ~ -^qPhadron + '^q ' (1^2) 

^G = ^GP,^ad.o„-AC^ (153) 

In terms of these variables, the hadron mass squared reads 

ML,^^ = ^^-t^^"^' + ^^^^ + MIg , (154) 

Xq XG 

where M'^^^^^j^g is now to be identified with the eigenvalue for the mass 
squared of the quark subsystem in a hadron and Alg is set to the corre- 
sponding eigenvalue Mq, which could be considered to have a ground state 
value or an excited value, if the condensed gluons were in an excited state. 
The constant is either 4m^ or 9m^ plus an appropriate number times 
muj with u = ojm or lo = ojb, respectively. Having introduced 

= T7^%y- . (155) 



Mq + MG' 

Mg 

Mq + MG' 



= T7^%r ' (156) 



one can define 




(157) 



= \ l^OXq - MqXG] . (158) 

y XqXG 
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Proceeding as in the case of two quarks in a meson, Eq. (112), one estab- 
lishes 



M 



qG 



+ 



(160) 



where two new parameters, (ph and M^, are introduced. 

The parameter ip^ is used in the similar way to how 'Pgiue is used in the 
invariant mass of quarks in mesons. However, ^p^ serves merely the purpose 
of notation for the unknown quantity of binding between quarks and G 
whose value must be adjusted in the process of correcting the candidate for 
the first approximation. In RGPEP, ip^ can be expected to depend on A. 
Intuitive arguments are offered in Appendix |Dj This candidate for a first 
approximation will only be satisfying if it turns out in future calculations 
in RGPEP that some optimal value of iph can be adjusted as a function of 
X/Aqcd and the corresponding coupling constant, assuming that the light 
quark mass parameters do not matter at small A where the effective quark 
masses are on the order of Aqcd anyway and do not vary with A so rapidly 
that no average constituent picture can correspond to QCD. Assuming that 
iph does not vary from hadron to hadron, one can introduce 



MgMg 

M„ + Mg 



(161) 
(162) 



The parameter /i/j is the reduced mass for the two-body system quarks-G. 
It depends on a hadron as far as Mq depends on a hadron, assuming Mq can 
be considered universal in the first approximation. Thus, the mass-squared 
eigenvalues for the whole hadron takes the form 



Ml = (Mg + Mcf 



Ml + 



Mg 



(163) 



where Uh denotes the excitation quantum number in relative motion of 
quarks with respect to the gluon condensate in a hadron, being zero in a 
ground state. Thus the momentum width of relative motion o quarks with 
respect to G is order iph- 

The mass parameter M^ is introduced by fiat and brings in another 
considerable degree of ambiguity. Physical motivation for M^ is that for the 
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gluon condensation to occur in hadrons it must be favorable energetically. 
If one just added the free energy of G and a harmonic quarks-glue potential 
energy to the energy of quarks, the condensation of gluons would only add 
energy to the system of quarks. The effect of condensation should rather 
be opposite: the concept of condensation of gluons inside hadrons implies 
thinking that the mass of a hadron state is lowered by condensation of 
gluons in comparison with a state in which condensation is absent. Since 
the kinematical minimum of energy of quarks and glue is Mq + Mq, which 
corresponds to = and Xq = f3g, the interaction mass parameter 
may be estimated by inspecting the condition 

m2 > {Mg + MGf + ^{Mq + MG)LOh-M^. (164) 

Assuming that the resulting lower bound on M^^ provides an estimate of its 
likely magnitude in QCD, we obtain 



Mi^^ M,Ma[2 + ^ + ^J, (165) 

where 3 signifies the number of spatial dimensions. Assuming that Mq equals 
nucleon mass mjv, and Mq ~ ~ uin/S, one obtains ~ 1.2 rn-Tv- If 
Mq ~ Mq ~ mjsf and Uh ~ w^v, we have M^p ~ 2AmM- This means that 
M^p can be adjusted to obtain a hadron mass as coming mainly from the 
quark mass eigenvalue Mq provided that M^ is quite sizable. Such sizable 
energy benefit from gluon condensation eliminates a large contribution of G 
to a hadron mass and sustains the possibility that a simple oscillator quark 
model without any glue can reproduce masses of light hadrons assuming 
that the effective quark masses are on the order of Aqcd 
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We wish to stress that Eq. (159) does not imply one and the same 



interaction between quarks and glue G in all hadrons even if the parameters 
concerning G are assumed the same for all hadrons. Different hadrons are 
actually having different Hamiltonians characterized by different quantum 
numbers of their quark component, such as radial and orbital excitations 
(elementary spin effects are discussed in the next Section) and the mass 
eigenvalue for the quarks, Mg, depends on these quantum numbers. This 
mass enters the definition of kh and thus also the definition of harmonic 
potential between the quarks and G that is defined in terms of d/dkh. As 
a result, potentials between quarks and G vary from hadron to hadron in a 



^* If the condensation mass advantage, —M^, were associated also with condensation 
of quark-anti-quark pairs, one might expect a pronounced reduction in masses of 
the hadron states in which a bilinear effective quark field expectation value plays a 
significant role in the dynamics, with tt mesons being the primary candidates. 
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well-defined pattern: the higher excitation of quarks, the greater their mass 
and the more momentum of a hadron carried by quarksj^ For example, 
quarks in excited nucleons are expected to carry a larger fraction of the 
resonance momentum than quarks carry in a ground-state nucleon. And vice 
versa, an excitation of the gluons condensed in a hadron increases Mq and 
the fraction of a hadron momentum they thus carry. Nevertheless, a CQM 
for light hadrons should be viewed as corresponding to n/^ = when A > Ac is 
actually lowered to Ac- As suggested in Appendix [D| when A is lowered to Ac, 
the role of G may be imagined taken over entirely by the content of extended 
constituent quarks that overlap each other heavily, so that no separate glue 
component can exist in hadron besides what is already contained in the 
extended effective quarks. Instead of the speculation, however, the proper 
problem is what will result from attempts to solve the RGPEP equation for 
the effective Hamiltonian at A = Ac. Strictly speaking, nothing is known 
currently about the solution. 



6. RGPEP in QCD and phenomenology 

The discussion that follows is limited to the case of all effective quarks 
having the same mass m. The common mass is expected to be a reasonable 
approximation for quarks u and d at A ~ Ac- Calculational complications 
due to differences in mass between these quarks are not discussed. 

Assuming that ^giue is on the order of Ag^-,^, which means that it has a 
similar value to the values obtained in QCD sum rules for '-p vacuums one ob- 



tains phenomenologically attractive values oIujm and ojb [HJ. Eqs. (112), 



(127), (154), (159), and (160), imply together eigenvalues and wave func- 
tions for light hadrons in the form of harmonic oscillator solutions that 
include the glue component G. The generic oscillator eigenvalue problem 
for two particles has the form 

io^ 1 \ 

— + -mu?r'^\i) = EiIj , (166) 
2m 2 J 

with eigenvalues En = {n + 'i/2)oj and the ground-state wave function tjjQ = 
A^exp [— p^/(2m-u;)] for n = 0, where p and fare canonically conjugated 
variables. 



For two effective constituent quarks in mesons, using Eq. (116) that 
explains what happens in the LF eigenvalue equation for A^gg in terms of 



a NR approximation to Aiqq, Eq. (115) implies the first LF approximation 
of the form 



M^-^ = 4m^ + 4(n + 3/2)mu;M , (167) 



Provided that the first-approximation quantity is kept constant. 
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V', 



99 



iVggOeXp 



1 



4ma;Af 



x(l 



4m 



(168) 
(169) 



with frequency um = T^^giue/ {^>f^) and wave functions of excited states, 
Tpqqm generated by building harmonic oscillator excitation operators from 
the vector k and gradient d/dk in a standard way and applying them to the 
ground state wave function ipqqQ. 



For three effective constituent quarks in baryons, using Eq. (145) that 



explains what happens in the LF eigenvalue equation for A4^g in terms of 



a NR approximation to Aisq, Eq. (144) implies the first LF approximation 
of the form 



+ 6(ni + 77-2 + 'i)muJB , 



iVsgoexp 



1 



+ m 



1 

Xi 



QmcoB 
1 

GmujB 

1 1 

— + — 

X2 X3 



9 



(1-X3)fc^ 
X1X2 



9 



+ 



X3(l - X3) 



(170) 
(171) 

(172) 



with baryon oscillator frequency ojb = \/5/8 ujm and wave functions of 
excited states, '4>3qnin2j generated in a standard way by building harmonic 
oscillator excitation operators from the vectors K and Q and gradients 
d/dK and d/dQ and applying them to the ground state wave function 
■03ijo- The LF momentum variables are defined in Eqs. (118), (120), (121) 
and ([T22J). 



In summary, the factor in a hadron ground-state wave function that 
depends on the relative motion of n quarks has the form 



'nqO 



NnqO exp 



1 



2nmuJr, 



[nm) 



(173) 



with n = 2 for mesons, n = 3 for baryons, 0:2 = com, ^3 = ^B, and pi the 
on-maas-shell four-momentum for quark number i. 

Exponentials of an invariant mass squared of quarks are popular as wave 
functions in phenomenological studies. Here such exponentials are related 
to gauge symmetry and reinterpretation of the gluon condensate as a part 
of a hadron. 



The spectrum of light hadron masses in Eq. (163), i.e., the spectrum 



corresponding to the ground state of relative motion of quarks with respect 
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to the glue G, becomes equal to the spectrum of masses of the quark com- 



ponent alone when the estimate ( 165 ) is adopted for the benefit of the gluon 



condensation in a hadron. This result reproduces success of the constituent 
quark models 



26 



At the same time, when the relative motion of quarks with respect to the 



glue part is described by Eq. (160), the hadron mass is given by Eq. (163) 



and the corresponding wave function is generated from the ground-state 
wave function 



1pqG0 



(174) 
(175) 



2{Mg + MG)iOh 

in a standard way for harmonic oscillators, first for the quark component 
and then, using the quark component mass eigenvalue, for the whole hadron 
if the quarks are excited in their motion with respect to G. Thus, the in- 
variant mass squared of quarks and G together, (Pg + Pg)"^ is evaluated 
in a standard way using their LF momenta and their minus components 
calculated using quark mass eigenvalue Mg and the glue mass Mc, respec- 
tively. Normalization factors are fixed by normalizing probability to 1, or 
fixing the charge to the appropriate value (if the charge is not zero these 
normalization conditions are the same). 



6. 1 . Form factors 

The relative motion of quarks with respect to the glue G smears the 
quark observables that follow from the quark factor in the hadron wave 
function alone. Consider the ground state of quarks' motion with respect 
to G. Calculation of the baryon form factor at a small momentum transfer 
q, = —Q^, < < Aq(-,^, is represented graphically in Fig. [ij Meson 
form factor calculation involves 2 instead of 3 quarks but otherwise proceeds 
in the same way. 

Using the Breit frame with g"*" = 0, which makes q~ that depends on 
masses irrelevant, one has = q^"^ and 

167r'^Xg(l — Xg) 

The oscillator functions of relative quark motion are independent of quark spins. A 
Coulomb part of the wave function must depend on spins. For small A, effective 
quark spin wave function can be treated as a separate factor in the oscillator part of 
the wave function. Spin factors can be constructed in a boost-invariant way using 
LF spinors and following the example of heavy quarkonia |18| in the case of mesons 
or loffe currents [3D] in the case of baryons. However, in contrast to models not 
related to QCD, RGPEP provides a scheme for systematic study of spin-dependent 
corrections to the first oscillator approximation. 
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Figure 1. The RGPEP calculation of a baryon form factor at A ~ Ac- 

X '4'*qGQ[Xq,l^q + {'^-Xq)q^]F^q{Q^)'4jqG0{Xq,K,g), (176) 

with normalization condition 

i6Ar^l*,c„K,s)^ (177) 



Using the wave function in Eq. (174) and changing variables to x = Xq and 
K,-^ = + {1 — x)q-^/2, one obtains 

FhiQ^) = Fsq{Q^)f{Q), (178) 



2 



16.3^(1 -x) '^^^°^^'"^''^" ^ (179) 
After integration over k"*", one is left with 

f(Q) = /Ml M^x ' (180) 



/o '^a;e 

where = (Mq + MG)coh- It is now visible that f{Q'^) is hardly differ- 
ent from 1 for small Q^: the inclusion of the glue component G does not 
significantly alter the result for a form factor calculated as if the glue com- 
ponent was absent. Moreover, in the picture discussed in Appendix [D| the 
glue component may disappear in favor of overlapping constituent quarks 
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Figure 2. The RGPEP calculation of a baryon structure function. 

when A is lowered below Ac- The momentum fraction carried by the quarks, 
X = Xg, becomes very close to 1 and the factor f{Q'^) becomes 1 

At large Q, much greater than a hadron mass, the hadron state needs 
to be represented in terms of the Fock components created by quark and 
gluon operators at A comparable with Q itself, in order to obtain a simple 
picture based on the smallness of an asymptotically small coupling constant 
that determines the strength of interactions which are responsible for dis- 
tributing the large momentum transfer g to a minimal set of constituents 
required to build a hadron. The minimal quark component may carry prac- 
tically the whole momentum of a hadron, Xq — )• 1. When this configuration 
dominates the transition amplitude, the whole hadron is turned from the 
total momentum P to P' = P + q in a combination of two mechanisms. One 
is — 7- 1, and the other is the distribution of q in the quark component. 
The former contributes to soft effects, and the latter is responsible for hard 
exclusive processes, cf. |4l|. 




The transition amp litude for deep inelastic lepton-hadron scattering is 
illustrated in Fig. |2P^| A hard photon or other boson is suddenly absorbed 




Even if the component G contained quark-anti-quark pairs, and were used to account 
for the quark condensate inside hadrons f20], its neutrahty would imply a small 
contribution to form factors of charged hadrons, such as proton. It might, however, 
contribute a detectable piece for chargeless hadrons, such as neutron. 
The lepton is not shown. Also, Fig. [Sjignores the possibility that the impinging boson 
is absorbed by a constituent inside the component G, which may be thought here to 
be only made of gluons. If G included quark-anti-quark pairs, it would contribute to 
the sea parton distributions as well as the quark condensate in hadrons, cf. [20] . 




6.2. Structure functions 
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by a constituent characterized by the scale Q. The new element RGPEP 
introduces is the possibility of using the transformation Wq\ of Section [4] to 
calculate the probability amplitude for finding such constituent in a hadron. 
The scale A refers to the particle operators in terms of which the hadron 
wave function is obtained from the eigenvalue problem. The scale Q refers 
to the hard boson. The final state in Fig. [2] is made of many constituents 
produced by Wq\ . Fig. [2] does not show that Wq\ actually acts on all 
constituents. 

The final state in Fig. [2] is in fact a virtual state whose evolution factor 
into observable particles is assumed to amount to unity. The energy, or P~ 
of the constituents that appear in the final state in Fig. [2j is dominated 
by the constituent at scale Q that absorbs the boson. This situation cor- 
responds to the leading operator expansion terms, i.e., hand-bag diagrams 
that are more important than cat-ears diagrams. One can use the Hamil- 
tonian HQ{hQ) to account for P~ of the virtual state. This Hamiltonian is 
actually the same as H\{h\). Therefore, Hq{Q) may count the energy of 
particles in the final state including the interactions that are responsible for 
grouping constituents at scale Q into spectators at scale A. Thus, the in- 
clusive sum over final states may be replaced by summing over a relatively 
small number of spectators at scale A and a potentially large number of 
constituents at scales between A and Q. The operator Wqx in RGPEP is 
hence expected to describe the evolution of structure functions with . 

The above qualitative reasoning must be verified by new type of calcula- 
tions using RGPEP evolution equations in place of other evolution equations 
in |54l [55t I56j . Since it is known that RGPEP incorporates the well- 
known splitting functions in QCD [25], there is no obvious reason for ex- 
pecting that RGEPEP evolution will significantly differ from known results 
where they apply. On the other hand, RGPEP provides the framework 
for combining the perturbative evolution with a non-perturbative hadron 
wave functions obtained by solving eigenvalue problem for Hx{b\). The 
eigenvalue problem describes saturation at small A. But the evolution pa- 
rameter A determines dependence of wave functions on invariant masses of 
constituent states. The invariant masses depend on transverse momenta 
and fractions x of longitudinal momentum in a specific way. One can thus 
expect that the evolution in x \b7\ [SH] is uniquely related in RGPEP with 
the evolution in Q^. 

6. 3. Other processes 

The approximate picture of a hadron at scale A > Ac as built from the 
effective quark component and from G treated as a scalar boson suggests 
thinking that the dynamics of G can be further approximated by suitable 
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effective Hamiltonian interaction terms. These terms would result in pro- 
cesses in which quanta of type G participate in strong interactions, softened 
by RGPEP form factors in vertices. Nothing can be said at this point about 
physical relevance of the approximate concept of quanta of the type G be- 
ing exchanged between hadrons. In particular, it is not excluded that an 
exchange of such quanta may contribute to diffractive processes in hadron- 
hadron collisions [59], or even in photon-hadron scattering [60] if a photon 
were allowed to contain its own glue component coupled to quarks. 



6.4- Connection with AdS/QCD 

The NR form of minimal coupling in LF Hamiltonians at small A appar- 
ently leads to variables and harmonic oscillator potentials that are similar 
to the ones that Brodsky and Teramond |49t I^Ul EI] identified as providing 
a correspondence between formulae for hadronic form factors in LF QCD 
and in field theory i n 5-d imen siona l AdS space [61^ l62] . The new var iable s 
k for mesons, Eqs. (106) and (107), and K and Q for baryons, Eqs. (123), 



(124), (125), (126), differ from standard LF relative t hree- momenta of con- 
stituents in the CRF For example, kgtandard in Eq. (108) differs in length 
and direction from k. But there is no difference between these momentum 
variables when they are small in comparison to masses. Thus, the new vari- 
ables are identical to standard variables in the NR dynamics at A > Ac. In 
addition, they fully describe a relativistic relative motion of constituents, 
always providing an exact representation of their free invariant mass, but 
in a different way than the standard momentum variables do. Standard 
variables kl 



^standard'' standard'' "^^^ Q standard 

ative momenta. The new variables k^ , 
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depend on the transverse rel- 
and , are defined using the 
constituent masses and fractions x of their total P"*", independently of the 
transverse relative momenta. At the same time, the transverse relative mo- 
menta are scaled by mass ratios and ■\/x{l — x). The square-root factor is 
precisely the one that appears in Brodsky- Teramond holography [ST]. 

To illustrate the possibility of correspondence between the reinterpreted 
gluon-condensate induced harmonic potential in LF QCD and a soft-wall 
potential in Brodsky- Teramond AdS/QCD holography, we use example of 
a meson form factor. The discussion is far from complete and raises many 
questions regarding interpretation of simple equations. However, it is needed 
for showing similarities and differences among different approaches. 

Let the struck quark and the spectator quark carry nearly whole hadron 
momentum. The glue component G provides additional smearing in the 



In LF dynamics, CRF differs from the bound state rest frame because conservation 
of P""" implies that is not conserved by interaction terms in LF Hamihonians. 
^° E.g., see Ref. [30]. 
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form factor expression, Eq. (176), unless it is totally absorbed in the con- 
stituent quarks Let us assume here that the latter option holds and all 
one needs to consider is two constituent quarks. 
In terms of the new variables, 



K^/^/4x{l-x), 

{x - 1/2) 2m/y^Ax{l-x), 



the LF eigenvalue equation for a meson of mass M reads 



(181) 
(182) 

(183) 



where = mcoM and r = id/dk. Using a factorized wave function, 
ip[k-^,k^) = N(l){k-^) fj^{x), where is an eigenfunction of the oscillator 
in z-direction, 



Hn{kz)e ^ 



H„ 



m{x - 1/2) 



K^yx{l — 

one arrives at the eigenvalue equation for (f){k-^), 
4m2 + 4:k-^^ + K^r^ ^ + (4r/ + 2)k^ 



i<P-x{'\. — x) 



(184) 



(185) 



which is the eigenvalue problem to compare with the Brodsky-Teramond 
holographic eigenvalue problem, e.g., Eq. (10) in [5T]. We denote quantities 
used by Brodsky and Teramond with subscript BT, except for their 



k' 



kQrp/2 . 



Eq. (185) can be written using kbt = V^k as 

.2^ 



- 4m2 - (2?? + , 



(186) 
(187) 



(188) 
(189) 



where k'^^ = — [d/dC,'^^'^ . For angular momentum around z-axis equal Iz, 

(190) 



using C = I C"*" I ) one has 

^2 1 d 



il 



See Appendix [d| 



48 SDGAPP2011GluonCondensate printed on July 1, 2011 



Defining (f. 



- 4'bt/VC, we arrive at 

1-4P 

^ + 4tC' + (2?? + 1)4t 



4(2 



(191) 



which by comparison with Eq. (11) in [5T] for massless quarks renders 

U{C) = ^^tC + (2r? + 1)4t , (192) 
as a counterpart of the potential C/(C) in Ref. [51], with 

27r 



4r 



2mujM 



(193) 



The role of quark masses on the right-hand side of Eq. (191) requires 
explanation. A sound explanation is currently not available in the sense 
that it is not clear why the quark mass is set to in Eq. (10) and (11) 
in Ref. [51J. But it is plausible for constituent quarks that the condensate 
mass-advantage constant in Eq. (163) can be considered a result of 
incorporation of G in the constituent quarks entirely. may be so large 
that it reduces the contribution of the quark masses to the hadron mass and 
thus reduces the term —Am'^ in the eigenvalue. 

On the other hand, the SW model eigenvalue Eq. (12) in Ref. [52] reads 



^P" + 
E 



4^ = Eip, 



in, + 2X, + 2. 



(194) 
(195) 



By dividing Eq. (191) by k^^-, and using variable z = kbtCi one obtains 



W{z) + 



li 



1/4' 



•}BT 



Ei 



'BT ■ 



E 



- - (2r? + 1)4t 



(196) 
(197) 



"BT 



Comparing this result with the SW model result, Eq. (194), and identifying 
Xz with Iz, one obtains 



W{z) 



(198) 



The SW model quadratic potential appears to have a coefficient given by 
the gluon condensate inside hadrons, Eq. (193). 
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If T/ = 0, which means no excitation along the LF, only —k^j- enters the 
eigenvalues. This is a sizable term (see below) and its compensation is not 
guaranteed by M^. 

Baryon form factors are described in terms of an active constituent and 
spectators. The analysis proceeds in a similar way to the case of a meson 
built from two constituents of different masses. This simplification is avail- 
able because the free LF invariant mass can be written in terms of the active 
constituent and the rest of a hadron as if the rest of a hadron was a single 
particle with its mass equal to the invariant mass of the rest of a hadron. 
The closest analogy to consider is that a baryon form factor and eigenvalue 
problem can be represented in terms of a quark and a di-quark. 

It should be stressed, however, that it is far too early for drawing a firm 
conclusion regarding connection between gluon condensation in hadrons and 
SW models on the basis of a pure RGPEP calculation. The reason can be 
seen using power counting [TB]. Besides uncertainties concerning x~ and 
inverse of id^ , spin-independent quark-anti-quark interaction term in a LF 
Hamiltonian density with a quadratic potential, which is absent in canonical 
QCD, may have a coefficient proportional to ^qcd^ compensate the 
dimension of x-*-^. There also exists a possibility [27] that a dimensionless 
square of a product r^P'^ is used instead of x^"^, where r is the relative 
position of arguments of two field operators on the LF, r+ = 0, and P is the 
momentum carried by a product of two quark fields, i.e., two constituents 
interacting by the potential. This possibility implies that the harmonic 
potential does not necessarily requires a factor Ag^-.^) because a gradient 

cancels dimension of a distance. Another Aq^^, may be needed to cancel 
the dimension of additional integration measure cPx-^. Still, the potential 
term must vanish when A — >• oo. If it vanishes as l/A'^, another Ag^-,^ 
may be expected in the coefficient. There always exists a possibility to use 
quark masses as dimensionful quantities. Together, one may need powers 
of Aqcd as high as 6. Such terms may be hard to establish quickly in any 
renormalization group procedure. 

In any case, the above results certainly allow one to entertain the possi- 
bility that the small-A RGPEP picture of hadrons, at strong coupling, can 
actually be the one that corresponds to the AdS picture, including the SW 
model |52| . On the other hand, even if the SW model does correspond to 
the Brodsky-Teramond LF holographic picture with a harmonic oscillator 
potential between an active constituent and spectators (in transverse direc- 
tions), and even if these pictures can be related to the gluon condensation 
in hadrons, how is RGPEP to explain emergence of the 5th dimension in 
the AdS/QCD analogy on the QCD side? 

It has been pointed out that the 5th dimension in AdS may correspond 
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to a renormalization group scale in gauge theory [631 l5^J • More recent cal- 
culations, such as [65, 66j, also point in this direction. We observe in this 
context that RGPEP provides an opportunity to explicitly introduce a 5th 



dimension in QCD with a simple physical interpretation. Namely, Eqs. (39) 



and (40) for effective fields can be rewritten using the observation that A 
corresponds to the inverse size of effective particles. The effective parti- 
cles exhibit this size in their interactions; the effective interactions contain 
vertex form factors whose momentum width is determined by the RGPEP 
parameter A. It is convenient to think about the size of the effective particles 
using variable s = 1/Ap^ 

So, instead of labeling the field and particle operators with A, we can 
label them with s that corresponds to the size of effective particles. The 
effective fields on the LF can be rewritten according to a rule 

^P,{x) = ^[k] bsik) e^^^^ = s) . (199) 

In this rule, the 5th argument of the quantum field is the RGPEP size 
of effective particles. This size plays the role of a renormalization group 
parameter in RGPEP. 

While the AdS dual picture involves its own renormalization issues [^7\ 
EHl ES ED], and RGPEP methodology may eventually find applications in 
these issues, too, the 5th dimension of particle size makes QCD a priori a 
5-dimensional theory, with scale invariance broken by Aqcd- The structure 
of such a theory may or may not be similar to an AdS field theory. On the 
other hand, RGPEP is certainly available for studying the 5-dimensional 
QCD. 

The issue that emerges immediately is that in order to evaluate observ- 
ables in QCD one in principle can use just one value of s or A. There is no 
need for integration over the 5th coordinate. So why would one integrate 
over the 5th dimension in the dual picture? 

The Brodsky-Terramond LF holography for hadronic form factors sug- 
gests that the integration over the 5-th dimension in AdS field theory cor- 
responds to integration over relative transverse momentum of an active 
hadronic constituent with respect to spectators. But RGPEP suggests that 



the LF wave functions have the form shown in Eq. ( 174), which is Gaussian 
in invariant mass of the hadronic constituents. Integration over the trans- 
verse relative momenta in form factors amounts to integration over invariant 
mass, with a proper rescaling by \/x{l — x). The invariant mass appears in 
ratio to a square of the width parameter A or, equivalently, in product with 



The particle size parameter s enters the non-perturbative RGPEP equat ions through 



replacement of the derivative d/dX on the left-hand side of Eq. (C.l I by d/ds 



„4 
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the square of the effective particle size parameter s = 1/A. The challenge is 
to explain how the integration over M can be turned into integration over 
A or s = 1/A. 

Although we do not provide here the ultimate response to this challenge, 
we offer a qualitative argument that suggests where one can look for the 
solution in future studies. Consider a fixed scale of an observable, such as 
Q in a form factor. For a given A = l/s, one can evaluate the observable 
in QCD by integrating wave functions over AA. But if the wave functions 
are functions of A^/A = sM, they are constant on hyperbolas defined by 
the condition that sM is fixed. One can consider a plane of variables M 
and s and plot a wave function in a perpendicular direction, forming a 3- 
dimensional plot of the wave function. Select a value of s = 1/Q and draw 
a profile of the 3-dimensional plot for the selected value of s as a function of 
Ai. Draw another profile as a function of s for Ai = Q. These two profiles 
are identical for all wave functions that only depend on the variables Ai and 
s through their product t. In such circumstances, integration over relative 
motion of quarks of fixed size Q is equivalent to integration over all sizes 
s of quarks that have fixed invariant mass Q. In QCD, this pure scaling 
picture for lightest quarks must be broken by the scale set through Aqcd- 
The small-A, or large-s picture of hadrons in RGPEP appears in ac- 
cordance with expectations that there exists a low-energy regime with a 

In the stabilization region, the gluon 
condensation parameter ipgiue inside hadrons is considered constant as a 
function of A. The hadronic wave functions are Gaussian in the invariant 
mass with the width provided by quark masses and ifgiue ■ The overall result 
is that at small A the effective parameters are frozen at values comparable 
with Aqcd- On the other hand, when A increases, and s decreases, the 
Hamiltonian changes and it is expected to eventually simplify to the QCD 
canonical form with counterterms at A — )• oo, i.e., for point-like quarks and 
gluons, s — )■ 0. Therefore, the wave functions of hadronic eigensolutions 
change with A, or s. One may expect that when the wave function width in 
momentum space increases and the sizes of effective quarks and gluons de- 
crease, the invariant mass dependence turns from Gaussian to a power-law 
behavior at large virtuality, perhaps as indicated by perturbation theory for 
exclusive processes in a coUinear approximation |41j . 

Finally, regarding the Brodsky-Teramond holography and SW models, 
we wish to mention that one considers different values of the SW parameter 
Av^ in front of in the AdS equations for masses of mesons (M) and baryons 



stable strong coupling constant |71 



RGPEP may be used to seek a precise definition of such coupling constant us- 
ing quark-gluon, three-gluon, and inter-quark potential terms in the corresponding 
Hamiltonians, cf. [381125) . 
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(B) [72]. Namely, km ~ 0.54 - 0.59 GeV and kb ~ 0.49 GeV, with 

— ~ 1.15 ±0.5. (200) 

KB 

In the reinterpreted gluon condensate image for hadrons, assuming a uni- 
versal value of mass for u and d constituent quarks, one has 

^ = f^V^'~ 1.125. (201) 

KB \5J 

This agrees quite well with the phenomenological result. 



7. Conclusion 

Reinterpretation of the gluon condensate in RGPEP has several impli- 
cations that matter in theory of hadrons. Instead of the entire vacuum, the 
gluons condense only inside hadrons. If this option were actually realized in 
QCD, there would be no need anymore to construct the quantum vacuum 
state in Minkowski space that satisfies fundamental requirements of invari- 
ance with respect to Poincare transformations, a construction that so far 
eluded theoreticians. At the same time, the LF Hamiltonian formulation 
of QCD would become free from any obligation to produce a non-trivial 
vacuum that is commonly assumed to provide physically important expec- 
tation values. When reinterpreted, the same expectation values come from 
distribution of matter that is limited to the interior of hadrons. Simulta- 
neously, it becomes not clear how one should interpret expectation values 
that are measured in lattice formulation of gauge theories; they might also 
correspond to the hadron interior rather than a state of infinite volume. In 
any case, the reinterpretation is available within the RGPEP that provides 
tools for a close inspection what actually happens. Thus, the picture of 
hadrons that is sketched here in the context of reinterpretation of gluon 
condensate can be viewed as providing an idea how the CQM may be an 
approximation to QCD and a method for verifying this idea. 

The constituent picture based on gauge symmetry dictates Gaussian LF 
wave functions that are exponentials of free invariant masses of constituents. 
The associated oscillator frequencies are in agreement with quark models 
phenomenology when the gluon condensate inside hadrons is set to the value 
found in phenomenology using QCD sum rules. The corresponding eigen- 
value problems appear naturally in terms of variables in which AdS/QCD 
models are developed, especially the SW model. The ratio of baryon and 
meson oscillator frequencies in SW models agrees with the ratio implied by 
the reinterpreted gluon condensate in RGPEP. 
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This article provides no comparable evidence for reinterpretation of the 
quark condensate. But one may hope that RGPEP can generate constituent 
quark masses when A is lowered to values comparable with Aqcd and shed 
this way some light on the mechanism of absence of chiral symmetry in 
hadronic spectrum and the nature of vr- mesons as nearly Goldstone particles. 

The author's opinion that RGPEP provides a method for verifying the 
reinterpretation of the gluon condensate and solving for hadronic struc- 
ture in QCD is mainly based on the fact that RGPEP is now available in 
both perturbative and non-perturbative versions, both being boost invari- 
ant. The step beyond perturbation theory maintaining boost invariance is 
seen as a chance for simultaneously taking advantage of knowledge about 
hadrons in the constituent picture and in the parton picture, while RG- 
PEP appears capable of generating required scale dependence dynamically. 
Seen this way, RGPEP provides an operator calculus for deriving quark 
and gluon wave functions of hadrons in the LF Fock space and unify cal- 
culations of the hadron mass spectrum and partonic distributions inside 
hadrons. At the same time, new variables identified in the reinterpretation 
of the gluon condensate, relative momentum three-vectors k in mesons and 
K and Q in baryons, provide a frame of reference in which LF solutions 
can be classified in terms of well-known angular momentum classification 
of constituent wave functions. The same transverse variables occur in LF 
AdS/QCD holography. Longitudinal variables are new and require further 
testing. 

A summary of RGPEP program for verifying the reinterpretation of 
gluon condensate is following. Start with canonical QCD. Set up regulated 
Hqcd on the LF. Apply RGPEP to lower the scale parameter A toward 
Aqcd- Find counterterms in Hqcd and evaluate effective H\. Solve eigen- 
value problem of H\ and see if the gluon-condensate induced oscillator pic- 
ture is reproduced by gluon components in hadronic states. So far, the 
RGPEP scheme is known to work only in some crudely approximate calcu- 
lations in QCD for heavy quarkonia. Studies of hadrons built in QCD from 
light quarks require breaking an entirely new ground. 

Appendix A 

Universality in RGPEP 

This Appendix explains the RGPEP scheme in which CQMs can be 
sought as limited representatives of the universality class which QCD is 
thought to belong to as a quantum field theory with a Hamiltonian acting 
in a Hilbert space, cf. [8]. The scheme is composed of two interrelated parts. 

In the first part of RGPEP, one constructs a sequence of rotations U\ 
for creation and annihilation operators in order to discover counterterms 
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in Hcxi and obtain the effective Hamiitonian H\. The second part invoives 
solving the eigenvalue problem of Hx. Solutions to the eigenvalue problem 
of Hx are needed to fix finite parts of counterterms in Hoo [73l [Tlj . 

The first part of the RGPEP calculation involves dealing with ultraviolet 
divergences. Instead of standard renormalization group procedure based on 
Gaussian elimination in linear problems [75^ I76j . RGPEP is based on a 
unitary transformation of field or particle operators. In terms of matrix 
elements of a Hamiitonian, this means that one rotates the basis states like 
in the similarity renormalization group approach (SRG) |73| I74j. 

Formulation of RGPEP in terms of operators for effective particles or- 
ganizes otherwise un-intelligible variety of Hamiitonian matrix elements in 
a physically motivated way; one traces coefficients of specific operators in- 
stead of only calculating matrix elements to which many operators may 
contribute. It also provides the concept of effective quantum field; see Eqs. 



(39), (40), and below. The operator structure of RGPEP results in uni- 
tarity features in the calculated interactions that are absent in the general 
formulation of SRG, since SRG applies also to Hamiltonians that cannot 
be written in terms of creation and annihilation operators. For example, 
Ux{boo) = Ux{bx)- The operator structure is also helpful in demonstrat- 
ing a connection between the momentum space formulation of the theory 
in terms of quanta (effective particles) and position space formulation in 
terms of quantum fields (effective fields). The interaction terms in effective 
Hamiitonian densities are non-local. RGPEP provides a method for calcu- 
lating the effective non-local interactions that correspond to renormalized 
canonical theories |27j . 

In contrast to the first part of the RGPEP scheme that is akin to the 
SRG scheme, the second part resembles the standard Wilsonian procedure 
because it involves Gaussian elimination, in the mathematical sense as part 
of solving a linear problem. However, in sharp distinction from the stan- 
dard procedure, the second step of RGPEP does not involve ultraviolet 
divergences. These are eliminated in the first step. 

The RGPEP differential equation for Hx in its generic form reads 

^nx = Txinx] , (A.i) 

where Tlx denotes Hx{boo) and J^x is a suitable functional '^^ The solution 

E.g., see [73- Appendix [c| provides a non-perturbative definition of TxlHx]. One 
needs to adjust tfie dimension of parameter A to dimensions of the Hamiitonian "Hx 
and functional J^x [Hx]. The parameter A is defined in terms of a momentum scale 
that plays the role of width in momentum-space vertex form factors, denoted also by 
A, or by the inverse of the width. The inverse of A corresponds to the size of effective 
particles that they exhibit in the interactions. The inverse is denoted by letter s that 
refers to the word size. 
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is thus also of the generic form 



fA 

■Hx = noo+ I dx':Fy[nx'] ■ (a.2) 



Hx{bx) is obtained by replacing boo in T~ix by bx- 

Irrespective of the theory one considers, for A much greater than the reg- 
ularization parameter A^/e (see Section[2]) times the largest allowed number 
of creation or annihilation operators in any product of them in the initial 
interaction Hamiltonian, J^x[T~(-x] is exponentially close to zero. The reason 
is that the regulated interactions are exponentially close to zero when the 
invariant mass of interacting particles is grater than allowed by the regular- 
ization. 

So, one can always introduce some A S> A/y^ and arrange a series of 
cutoffs Ao = A, Ai = A/2, A2 = A/4, etc., until one reaches A = A„ = A/2" 



for some large n This is done in analogy to Refs. fTS^ i76| for the purpose 
of sequencing the RGPEP procedure of evaluating Hx into manageable steps 
that are free from huge terms resulting from ultraviolet divergences of the 
local theory. Namely, step number k in the sequence that builds up the 



integral in Eq. (A.2), contains only integration over scales between A^+i to 
Afc. When one evaluates effects of interactions that allow changes of invariant 
masses only in the range between A^+i and Afc, the ratio Ao/An — )• 00 does 
not contribute. Moreover, one can carry out each and every one of the finite 
steps using well-known techniques, such as perturbation theory or numerical 
methods. 

Although the sequence one obtains is analogous to the sequence consid- 
ered in Wilsonian renormalization group procedure [76], it differs consider- 
ably because no states are eliminated from the domain of the Hamiltonian. 
The RGPEP steps are automatically free from the dependence on eigenval- 
ues that appear in Gaussian elimination in matrix notation for eigenvalue 
problems and there is no need for an additional (and in principle deter- 
mined only up to a nearly arbitrary rotation of basis) procedure required 



Instead of the factor 2, one can choose e or any other similar number. The exponen- 
tial spacing of As in the sequence of cutoffs is introduced for the purpose of handling 
logarithmic singularities. In the case of simple models in which one can easily intro- 
duce an exponential grid in momentum space, such as the NR Schrodinger equation 
for one particle in a (5-function potential, it is most convenient to use a grid in which 
the chosen factor for cutoff reduction, similar to 2 or e above, is also equal to an 
integer power of the constant used in the grid, i.e., if = poa'' in the grid, then 
An/A„+i = a' with an integer I such that a' equals 2, e, or some other convenient 
factor chosen for the cutoff reduction. This relates the value of a chosen for the 
exponential grid in momentum variables with the RGPEP cutoff reduction factor in 
terms of an integer I. 
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for maintaining formal hermiticity of effective Hamiltonians 

The ratio of subsequent parameters, A^+i/Afc does not depend on k and 
one can seek regularity in how successive integrations transform Hx^. to 
Hx^^^. Every element in the resulting sequence of Hamiltonians can be 
written in terms of dimensionless momentum variables = p'^/Xk instead 
of momenta p"*". Creation and annihilation operators also require rescaling 
in order to keep their commutation relations in terms of variables inde- 
pendent of the step number k. In principle, such rescaling corresponds to 
scale-dependent renormalization constants for fields in standard perturba- 
tive formulation of quantum field theory. The coefficients cx^ of products 
of creation and annihilation operators bx^, in Hx,, can be studied as func- 
tions of 2/"*" and x. These functions evolve in k. Their evolution can be 
classified in terms of characteristic dominant behavior associated with uni- 
versality classes, presumably associated with fixed points, limit cycles, and 
even chaotic behavior. Such analysis is also expected to help in identify- 
ing dominant effects due to small-x cutoff parameter 6 in gauge theories 
quantized on the LF (see Section [2]). No farther comments are offered here 
regarding the evaluation of Hx, except for stressing that the establishment 
of finite parts of ultraviolet counterterms eventually requires a reference to 
the second part of the RGPEP scheme. 

The second part involves solving for the spectrum of Hx„- This is fa- 
cilitated using techniques similar to the ones used in Ref. [78] with the 
following qualitative distinction. The Hamiltonian Hx„ = Hqx„ + Hjx„ is so 
narrow on energy (invariant mass) scale (defined by eigenvalues of Hqx„), 
that the wave functions of its eigenstates are typically exponentially sup- 
pressed outside the range of energies (invariant masses) that are similar in 
size to the observable energies (invariant masses) one is interested in |79j . 
Therefore, one is no longer in need to solve the severe ultraviolet renor- 
malization problem when carrying out calculation that removes residual 
(exponentially small) cutoff effects using a procedure such as in [78] . This 
suggests, for example, that the non-linear operators that produce correc- 
tions to scaling in Wegner's sense [80] only lead to small corrections to the 
scaling properties (if any are obtained) that result already from the first 
step of the RGPEP scheme. 

The expectation of lack of significant sensitivity of the second part of 
RGPEP scheme to the actual values of smallest eigenvalues can be forecast 
by analogy with (and it is confirmed in) an elementary model of a harmonic 
oscillator with an additional potential proportional to the fourth power of 
distance (a Higgs-like field dynamics in zero dimensions) |81j. Excited states 
of the pure oscillator (excited by multiples of hco) serve as a model of the 



In Ref. [7^, the operation chosen to maintain Hermiticity is denoted by R, Eq. 
(4.1). 
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Fock space basis with different numbers of effective particles, each of mass 
m = huj. The Hamiltonian with the quartic term models Hx^ after the 
entire first part of RGPEP scheme is completed. Namely, the interaction 
is narrow on the energy scale because it is able to change the number of 
particles only by 0, 2, or 4, which means that the size of energy changes due 
to interactions is limited from above by 4:huj, or 4m. The result of numerical 
computation of eigenstates in that model [81j is that for obtaining smallest 
eigenvalues accurately it is sufficient to use Gaussian elimination assuming 
that the smallest eigenvalues are simply zero. One can come down this way 
to cutoffs on energy of basis states not much greater than about 10 times 
huj, or 10m. This holds even for quite large values of the coupling constant 
that characterizes the size of the quartic coupling. 

Since the entire RGPEP procedure is in fact a sequence of successive 
approximations with an increasing space of variables when regularization is 
lifted, both parts of it, the rotation and solution, need to be iterated, in 
the sense of Wilsonian triangle of renormalization, until a stable result is 
established. The iteration involves finite parts of counterterms. They are 
fixed by comparison with data in the solution part. But the finite parts 
influence the rotation in the first part. This is why the two parts of the 
procedure are interrelated and neither the part using unitary rotation of 
particle operators nor the solution part with elimination of states can be 
distinguished as determining the path of successive approximations for the 
Hamiltonian on the triangle. As a result, the evolution of a Hamiltonian 
on the triangle is neither purely unitary nor a plain solving of an a priori 
specified eigenvalue problem. In other words, the ultimate self-consistent 
solution to the triangle of renormalization defines the theory one is interested 
in solving to explain observables. 

RGPEP requires extensive studies in order to verify if it can produce 
a universal low-energy theory that explains success of the CQMs starting 
from the Lagrangian for QCD. In addition, all of these models could also be 
subjected to RGPEP or SRG as long as they can be written in a Hamiltonian 



form, Therefore, there arises a question if all CQMs that do reproduce 
a considerable number of observables for smallest-mass states in hadronic 
spectrum actually converge on a single, universal model with some specific 
shape of quark potential. Then, the question would be if a universal model, 
if it is obtained, matches the effective theory obtained from QCD using 
RGPEP. 

It should be stressed, however, that the procedure outlined here is not 
limited to RGPEP for quarks and gluons in LF QCD. It can be applied to 
studies of universality in all Hamiltonians that can be expressed in terms 



This is not immediately obvious in the case of approaches based on the Dyson- 
Schwinger and Bethe-Salpeter equations. 
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of creation and annihilation operators acting in some Fock space and, even 
more generally, to all Hamiltonians that require renormalization and can 
be subject to the SRG procedure [73l [71]. In particular, one can apply the 
same procedure when using beautiful Wegner's flow equation [82} [83 1 IMII85] . 
and the latter can also be suitably altered in order to improve weak-coupling 
expansion that may apply in the case of asymptotically free theories [iSl |l6] . 



Appendix B 

Calculations of W 



Assuming that Aqcd formally tends to in comparison to quark masses, 
so that gx can be treated as extremely small when A is comparable with the 
masses, W could be calculated in RGPEP using perturbative formulae given 



in Section 4.3 and expanding LF quantum fields into bare creation and an- 
nihilation operators. In this case, the dominant term besides 1 is gxWi. 
However, for a realistic value of Aqcd, which is much larger than u and d 
quark mass parameters in the SM, and for A ~ Ac that is expected to corre- 
spond to the CQMs, perturbative expression cannot be considered reliable 
on the basis of smallness of gx. Nevertheless, perturbative calculations de- 
scribed here illustrate the structure of W that carries over to values of A 
comparable with Ac for as long as Tioo as an initial condition in RGPEP 
equations is replaced by the corresponding Hx and A is sufficiently near 
in magnitude to Ac so that W is not violently different from 1. The key 
argument is that the strength of the interaction is not fully determined by 
gx alone. Namely, when the RGPEP form factor fx is narrow as function 
of momentum variables, the net strength of the interaction is suppressed 
and the interaction can be weak even if gx is sizable. Thus, the calculation 
given below illustrates not only how the lowest-order perturbative W looks 
like in the case of small gx but also how the structure of W is related to 
the structure of T-lx that is weak because A is small, even though Tlx is 
expected to significantly differ from Tioo when A approaches Ac ~ Aqcd for 
light quarks. 

For the canonical quark field, we have on the LF 

= Yl /[^] [xcUfkabkacfe''^'' + XcVfkadl^^fe''''''j^^^^ . (B.l) 

acf 

The momentum and spinor notation involves (e.g., see Ref. |18j) 

[k] = e{k+)dk+d^k^/{16Tr^k+), (B.2) 
Ufka = B{k,mf)uf^, (B.3) 
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Vfka = B{k,mf)vfa- 
Vf„ = Cu}^, 

Ufa = V^2m7 C/ Xa 



(B.4) 
(B.5) 
(B.6) 



where C/ is a normalized vector in flavor space, Xa is a spinor that corre- 
sponds to one of two spin states of a fermion at rest, both normalized to 1 
and orthogonal to each other, C is the charge conjugation matrix. 



B{k, m) 



Vk+ 



m 



-.[A+k+ + A^{m + k^a^)] 



(B.7) 



is the LF boost matrix, and A± = 2 7o7 • For the canonical gluon field, we 
have 



'fco- 



(TC 



c^M „. „-ikx . fJ,* \ ikx 



ka ~ 0' ^fccr 



2k^ei/k-^,ej^). 



x+=0 



(B.8) 
(B.9) 



The same expansions can be used for effective quantum fields in which 
creation and annihilation operators correspond to some scale A [2ZF^ see 



Eqs. (39) and (40). 



= Hqoo- According to Eq. (32), 



To evaluate Wi using Eq. (31), one first needs to evaluate Hiooi and 



Ooo 



'Hiool = dx d X hiooi ■ 



(B.IO) 
(B.ll) 



Using Eq. (33) and the field expansions, one obtains the normal-ordered 



expression for Hq, 



Tin 



E/w 



(B.12) 



In order to relate the above notation to the one introduced in Ref. [27] , one needs to 
treat with positive fc^ as b with negative fc^, and with positive as a with 
negative k'^ , taking into account the well-known constraint relations for fermion field 
V^- = A-V- and gauge boson fields A' gH |86l |87l IM US] . 
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Similarly, using Eq. (34), one obtains Hi: 



ool ) 



123 



+ ^123 44^3 + h.C. 



(B.13) 



where Pc denotes the total momentum of created particles, Pa denotes the 
total momentum of annihilated particles, rA5(l,2,3) denotes the regular- 
ization factors described in Section [ij t?- = xL^^Xjc 



123 



if 



C1C2C3 



£162 • - ei£3 ■ e2K- 



2:2/3 



£263 • ^iK- 



^1/3 



,(B.14) 



with £ = and k = k^ — k^k^/k^, and h.c. denotes Hermitian conjugation 
of the 4 explicitly written terms in the square bracket in Eq. (B.13). 

The result for T-Ljooi needed in Eq. (31 ), can be written in an abbreviated 
form 

niooi = Yl [[i^23]2{2Trf 6\Pc-Pa)r Asil, 2,3) {1,2,3} , (B.15) 

123 

where {1,2,3} denotes the 4 + h.c. terms with 600, doo, and Ooo replaced 
by ^A, dx, and ax, respectively. Using this notation, the exact perturbative 
result for Wi in RGPEP reads 

Wi = Y J[l23]2{27Tf6^{P,-Pa)rA5{h2,3) 

{1,2,3}, (B.16) 



123 

fx - fx. 

Pa - Pc 



where and P^ are eigenvalues of Hq and / denotes the RGPEP form 
factor, 



fx 



(B.17) 



in which the squares of invariant masses of created particles, M."^ = P^, 
and annihilated particles, ^A'^ = P^, are evaluated using eigenvalues of T-Lq 
associated with the particle momentum components kf and A;^ for i = 1,2,3 
in every term [27] . 

The structure of {1,2,3} implies that Wi can replace a quark by two 
particles: a quark and a gluon. It can replace a gluon by a quark-anti-quark 
pair, or by two gluons. Reversed changes are also possible. 
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According to Eq. (30), evaluation of W2 involves action of Hiooi twice 
and T-Lioo2 once. This means that W2 can replace one virtual particle by 
3, or vice versa, change the number of virtual particles by 1, or do not 
change their number at all but only alter their individual momenta and 
other quantum numbers as dictated by Tiioo- 

When one evaluates W for A and Ac using expansion in powers of gx, 
creation and annihilation operators in Eq. (B.16) being those corresponding 
to A, the coefficient ifx — fxJ/{Pa—P^) only appears with first power oi gx. 
Higher powers involve different coefficients and the entire sum of the series 
in powers of gx is not known. Moreover, a complete answer may contain 
dependence on Aqqd that perturbation theory cannot identify. Therefore, 
the question is what structure of W one can expect for A ~ Ac ~ ^qcd- 
In particular, seeking a connection between CQMs and QCD, one needs to 
estimate the structure of W that is responsible for changing the number of 
effective particles, since CQMs do not include such interactions. 

When A is near Ac, one knows that W is near 1 irrespective of the size 
of both As and gx, because is 1 when A = Ac- However, the deviation of 
W from 1 is not the same for A ~ Ac ~ ^qcd as in the perturbative case 
discussed above for two As that are much greater than Aqcd- The significant 
change in W in comparison to the perturbative result for large As comes 
from the difference between Tioo and full solution for Tlx- Nevertheless, if 
one assumes certain particle-number changing interaction in Tix for A near 
Ac, W will put this structure of Tlx to action in a similar way to the one 
described above. The reason is not the smallness of gx but the presence 
of the form factors fx- They weaken the interaction by limiting the range 
of momenta that interacting particles may have. Thus, one can look at a 
candidate for Hx for some value of A near Ac as an initial condition for 
integration over a relatively small range of A near Ac and treat the whole 
interaction that changes the number of particles as weak even if gx appears 
large. 

The weakness of particle-number changing interactions with small A is 
reinforced by the particle masses as described in Section [2j Both quark, m, 
and gluon, nig, mass parameters in the region of A ~ Ac may be sizable. 
Gluon mass terms for effective gluons are not excluded by local gauge in- 
variance of Lagrangian for QCD because the interaction terms in LF QCD 
Hamiltonians with finite A are not local and their non-locality is not of 
canonical type. Since regularization of LF QCD Hamiltonian requires gluon 
mass counterterms and the latter have finite parts that are expected to de- 
pend on Aqcd, perturbation theory for scattering in femtouniverse [89] can- 
not be used to argue that effective gluons cannot have sizable mass terms 
for A ~ Ac. The spectrum of masses of lightest hadrons, which is quite 
sparse in comparison to the QED near-threshold atomic spectra with mass- 
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less photons, suggests instead that there is a considerable mass gap involved 
in excitation of gluon degrees of freedom. Since QCD is assumed to describe 
these data, it seems reasonable to expect that RGPEP equations in QCD 
lead to massive effective gluons, instead of massless ones. 

As an illustration of what one might expect of W in QCD in the range 
of A ~ Ac ~ ^QCD when one reasons by analogy with the perturbative solu- 
tion at large A, consider the structure of local three- gluon interaction term 
that originates from the gauge-invariant Lagrangian density —TrF^^ F^i,/2 
in canonical theory. The interaction causes splitting of a bare gluon into 
two bare gluons. Its structure is given in Eq. ( |B.13 ), as the term with Yi23- 



In distinction from this local canonical term (in which A = oo), the corre- 
sponding term in an effective Hamiltonian at small A ~ ^qcd, is non-local 
|27j . Suppose the complete non-local three-gluon vertex includes the same 
structure 1123 and a vertex factor Va(1, 2, 3) times the form factor fx. Such 
structure is known to emerge through order gf^ [25]. The corresponding 
l-to-2-gluon term in W of lowest order in the effective interaction reads 



Wy121 = Yl [[m'2{^T^fSHPc-Pa)^_ 



fx - fx 



123 

-.t J 



Pc 

X yi23l6.(l,2,3) a{^a{^axz, (B.18) 

where Va(1,2,3) plays the role of initial condition, instead of ^00(1, 2, 3) 
that provided the initial condition in T-Loo- Note that now the regularization 
factor rA5(l,2,3) is absent. This follows since the form factors f\ and f\^ 
entirely eliminate the ultraviolet divergences and they eliminate the small-x 
divergences when one adopts the assumption that effective gluons are heavy 
|16j . In these circumstances, the regularization factors are immaterial. 

Suppose Ac is smaller than the gluon mass nig at Ac. In this case, for 
A '-^ Ac, one can approximate P~ by m^/pj^ and P~ by Am^g/p^ in the 
difference in denominator while the form factors can be approximated using 
(see [27| for details) 

/, = ^-[Mh-^^\'/>^\ (B.19) 

with M-12 — ^{1^1 + ^^)) rewritten as 

fc^ + 3m^/2 g2 

f^ = e-(3"*?A')'e w5j^(W. (B.20) 

For A smaller than rUg , the relative momentum | /c | of gluons 1 and 2 is 
smaller than A/2 < mg/2 and the form factor can be very well approximated 
by Gaussian 
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The exponential in front causes the form factor to vanish when A is much 
smaller that discussed in Section [2[ This exponential is a part of 

the mechanism by which the presence of RGPEP form factors suppresses 
interactions that change the number of effective particles. The other part 
of the mechanism is the Gaussian factor whose width decreases when A 
decreases. This factor causes the range of interaction in momentum space 
to decrease and thus to diminish the net effect of it. For the form factor 
suppression mechanism to work, the coupling constant g\ must not blow 
up to huge values and compensate for the smallness of the form factor. 
However, such compensation should not be expected in QCD if the theory 
is to explain the phenomenological success of the CQMs. 

For Ac much smaller than m^, one could literally set fx^ to in Eq. 



(B.18). The resulting estimate for the structure of Wyi21 reads 

24m2fc2/A4 



Wyi21 = e-'</^'Yl l[mm^fS\Pi2-P3) 



e 



3mj/pJ 



123 9 
^A1"A2 



X Yu3Vx{l,2,3)a{.a{,ax3. (B.22) 



This term in Wi replaces a gluon with a pair of gluons. The relative motion 
of the two emerging gluons corresponds to a Gaussian wave function of the 
form e"'^*^^/'^^ with c = 24m^/A^. Assuming that A is comparable with 
nig ~ 1 GeV, the width of the Gaussian wave function is about A/5 ~ 200 
MeV, which is not unreasonable. Nevertheless, one has to remember that 
the full vertex contains also the function Va(1,2,3) which causes that the 
non-locality of Wyi21 is not determined solely by the Gaussian factor in 



Eq. ( |B.22D . 

Another illustration is provided by consideration of the first term in 
{1,2,3}, i.e., the term that replaces a quark of scale A by a quark and a 
gluon of the same scale. The only new elements in the analysis of this term, 
in comparison to the analysis described above in the case of gluons, are 
different values of masses and spinor factors instead of Y . However, the 
conclusion is similar: W changes a quark into a quark and a gluon. 

The larger the difference between A and Ac, the more important particle- 
changing interactions in H\. The full solution for W may replace a state 
of two or three constituent quarks by a superposition of states with many 
different numbers of gluons (and additional quark-anti-quark pairs created 
from additional gluons) at A, as dictated by the interaction terms in H\. 

The perturbative calculation of W , and calculations for small difference 
between A and Ac using assumed simple picture at Ac and guesses concerning 
interactions at A > Ac, can be replaced by non-perturbative solutions of 
RGPEP equations, see Appendix [Cj 
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Appendix C 

RGPEP beyond perturbation theory 

For the principles of renormalization group procedure for Hamiltonians 
|73j to apply in a boost-invariant LF dynamics, the right-hand side in Eq. 



(A.l) must guarantee that matrix elements of H\ vanish when changes of 
the invariant masses exceed A and become comparable with the invariant 
masses themselves. This feature will be called narrowness in invariant mass 
of interacting effective particles of width A, or just narrowness. One way of 
securing narrowness of width A in perturbation theory is to work with vertex 
form factors in RGPEP as described in Appendix [Bj While this option is 
welcome in perturbative calculations [l6], it involves a derivative of the 
Hamiltonian in the definition of functional Tx and the derivative is not easy 
to calculate outside perturbation theory. The perturbative approach of Ref. 
|74j defines differential equations in which the functional T\ guarantees 
narrowness of width A without depending on derivatives of the Hamiltonian 
on the right-hand side, but the generator of the transformation is again 
defined recursively in terms of a commutator of the generator with the 
Hamiltonian. The recursion formally secures an expansion to all orders 
in the coupling constants, but it is difficult to solve beyond perturbation 
theory. 

Non-perturbative formulation of RGPEP for operators in the Fock space 
can be defined in the form of a double commutator ^SOj [911 |92j equation that 
reads 

-^^x = mox^UpxlUx] . (C.l) 

The Hamiltonian Hx is obtained when boo in is replaced by bx- The 
required transformation Ux is thus meant to be always incorporated with 



accuracy dictated by the accuracy of solving Eq. (C.l). Strictly speaking, 
nothing is known yet about the accuracy one can actually achieve this way 
beyond perturbative RGPEP in QCD. 



The structure of Eq. (C.l) resembles the beautiful equation for Hamil- 
tonian matrices that Wegner proposed for solving theoretical problems in 
condensed matter physics [82l EH |8l]j^ In Wegner's proposal, is a 
diagonal part of the Hamiltonian matrix, 1-Lix is the off-diagonal part, and 
Tipx is the matrix Hix itself. Attempts have been made to apply the Weg- 
ner equation to Hamiltonian matrices obtained from quantum field theory 
for the purpose of eliminating matrix elements that involve a change in the 



To avoid a misunderstanding, it should be clarified that in his original work Wegner 
did not relate his equation to renormalization group procedures in condensed matter 
physics or quantum field theory. 
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number of virtual particles. Such matrix elements are found in the ini- 
tial conditions set using interaction terms from canonical Hamiltonians in 
quantum field theory. The Tamm-Dancoff truncation to a few Fock compo- 
nents was used to define a sufficiently small Hamiltonian matrix for using 
Wegner's equation [931IM1 l95l l96] . 

Wegner's equation with matrix 'Hp\ = Hix implies that the matrix Hx 
with small A is narrow with a width order A on energy scale (it would be P~ 
scale in LF dynamics). However, Wegner's equation for Hamiltonian matri- 
ces does not agree with kinematical LF symmetries, including the Lorentz 
boosts that are essential for explaining a connection between the constituent 
picture and the parton model picture of hadrons. Also, one desires an equa- 
tion for a Hamiltonian operator that in principle may act in the entire Fock 
space rather than only in a severely truncated space in a Tamm-Dancoff 
approach. These two issues will be addressed below by defining the oper- 
ator T-Lp\ and showing that the resulting equation leads to narrowness of 
1-L\ that respects all kinematical LF symmetries and is narrow in the entire 
Fock space. 



Appendix C.l Boost invariance 

The origin of the difficulty with LF boost symmetry in Wegner's equa- 
tion is that the left-hand side in his equation is supposed to be a derivative 
of a Hamiltonian with respect to a parameter, while the right-hand side is 
tri-linear in a Hamiltonian, an object clearly depending on a frame of refer- 
ence in a different way than a Hamiltonian does. LF power counting [IgJ for 
Hamiltonian terms suggests that a straightforward application of Wegner's 
equation in the SRG scheme may create a host of complex counterterms 
required for restoring kinematical LF symmetries. An attempt to cure the 
situation was undertaken [93 [98] by assuming a double-commutator equa- 
tion for matrix elements of the invariant mass squared rather than a Hamil- 
tonian. It was assumed that the LF boost invariance could be maintained 
because powers of the invariant mass are invariant with respect to boosts if 
the mass itself is. However, invariant masses of Fock states depend on spec- 
tators, resulting in effective interactions that depend on spectators. This 
means that the resulting effective interactions require extra care to recover 
the cluster decomposition principle [9.9| and it is not clear that they can help 
in literally solving quantum field theory where this principle is expected to 
be valid in effective theories, including effects of a complete renormalization 
procedure. 

The issue of LF symmetries in RGPEP was considered before [23]. Here, 
the operator Tipx is defined using the Hamiltonian Tix- Assuming that T-Lx 
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conserves momentum and is of the form 

= E E [1-^] ^^(1' ^) ll-'-dn, (C.2) 

where q denotes annihilation operators boo, the operator Tipx is defined to 
be of the form 



n=2~ V k=l J 



+ J 



ql - ■ -Qn, (C.3) 



and the symbol ^ denotes integration over momenta and summation over 
discrete quantum numbers that label creation and annihilation operators 
while p'^ denotes the +-component of momentum that labels the creation 
or annihilation operator number k. In words, Tipx differs from T-L\ by mul- 
tiplication of its terms by the square of +- component of total momentum 
carried by the particles that enter a term, which is the same as the +- 
component of total momentum carried by the particles that leave the term. 



In the resulting Eq. (C.l), both sides behave in the same way with re- 
spect to operations of kinematical LF symmetries, and the RGPEP width 
parameter A is invariant with respect to boosts of kinematical LF symmetry. 



The operator Hqx in Eq. (C.l ) is now set equal to a free-particle Hamil- 
tonian Tio, i.e., an operator built from products of one creation operator 
and one annihilation operator per particle species. The notation TIqx = Hq 



also indicates that Hqx in Eq. (C.l) is chosen here to be independent of 



A 40 It is a diagonal operator in the Fock space spanned by states created by 
operators 6oo from the LF vacuum state. Diagonal matrix elements of T-Lq 
are not equal to the diagonal matrix elements of the full Hamiltonian. Typ- 
ical Hamiltonians include interactions that contribute to diagonal matrix 



elements 



Appendix C.2 Non-perturbative narrowness in the Fock space 



Eq. (C.l ) is designed for operators that can act in the entire Fock space 



(action on every state is well-defined). Eq. (C.l ) with a constant Hq renders 
A-dependent interaction terms that die out when the change of the invariant 
mass squared (evaluated using the kinematical momenta and eigenvalues of 
Tio) of particles involved in an interaction exceeds A (spectators do not in- 
fluence Hamiltonian interaction terms). The narrowness feature is obtained 



This assumption simplifies the analysis that follows. 

It is possible to include interactions in Hox in Eq. (C.l I in the form of mass effects 
(self-interactions) and beyond (potentials), but these options are ignored here for 
simplicity. 
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through a universal mechanism for double commutator evolution equations. 
The factors of in T-Lpx do not change the general mechanism. 



We start from the observation that Eq. (C.l) can be solved using a 
method of successive approximations. An approximation is defined by how 
many terms are kept in the Hamiltonian. In principle, a Hamiltonian with 
finite A contains infinitely many terms. The terms that are kept in an 
approximation can be defined by the condition that they contain no more 
than a prescribed number N of creation and annihilation operators in a 
product. The number labels the approximation. Terms that contain 
more operators in a product are ignored in the approximation. Successive 
approximations are labeled by increasing numbers N . 

Suppose one wants % to contain at most N creators and at most N 



annihilators,^ Such Hamiltonian will be denoted by T-Lxn- Using notation 



corresponding to Eq. (C.2), this means that 



N ^ /■ 

H\N = X] ^^['^■■■^c + na]hxNn,na{'^,-,nc + na) tn.naA^-^) 
nc=l na=l 



where 



ric+na 



k=l l=nc+l 

Uc is the number of creation operators, Ua is the number of annihilation 
operators, and subscripts denote also all relevant quantum numbers of the 
operators they label. 

Suppose one wants to know the terms in T-Lxn that contribute to the 
dynamics of a physical state that certainly contains a specified Fock com- 
ponent IV') of the form 

. n 

|V> =|;[l...n]V'p(l,...,n) llqi\0), (C.6) 

k=l 

where P denotes the fixed total kinematical momentum of the state and 
otherwise the wave function ipp{l, ...,n) is not known. For example, in 
the case of mesons one may select the state that contains a quark and 
an anti-quark, and in the case of a baryon a state that contains three 
quarks. But one could also consider a state that contains 1182 quarks and 



''■^ One can limit the number of particle operators for different species of particles dif- 
ferently. In this case A'^ becomes a vector with natural components corresponding to 
species. 
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many gluons and quark-anti-quark pairs in order to describe a collision of 
two nuclei of gold neglecting all interactions but the strong. 

Having specified the set of operators that are included in Hxn, with 
ric < N, Ua < N, and unknown coefficients hx^ucUa^ can enumerate 
forms of all possible states that can be generated from by action of 
T-LxN on it once. All these states are eigenstates of the total kinematical 
momentum operator with one and the same eigenvalue P that characterizes 
Knowing the set of states that are obtained from by acting on it 
with HxN once, one can construct the set of all states that can be obtained 
by acting with HxN on twice, and so on. 

We introduce the set of states, denoted by R, that can be generated 
from I'i/^) by acting with HxN on it r times and are normalized in a limit of 
infinite volume. Then we introduce a minimal subspace in the Fock space 
(a space of smallest possible basis) that is sufficient to build the set R. This 
subspace in the Fock space is denoted by RripN- For example, if one starts 
with a state of 3 quarks, action just one time by the T-LxN that is equal to 
regulated canonical Hamiltonian of LF QCD, r = 1 and = 3, on |3g) 
produces states with 3 quarks, 3 quarks and a gluon, 3 quarks and a quark- 
anti-quark pair, and 3 quarks and 2 gluons. The space -Risgs is built from 
the Fock space basis states that are needed to construct states with these 
particles. Action twice, r = 2, generates additional gluons and quark-anti- 
quark pairs. The corresponding space R2 3q3 is built by adding the Fock 
space basis states that are needed to construct the additional states. When 
one assumes that TixN has more terms than the canonical LF Hamiltonian 
for QCD [l6], the space RrtpN is greater than in the canonical case. States 
created from j-^) with k particles by r actions of Hxn may contain up to 
k + t{N - 2) particles. 

The set R of states one obtains for any r is a priori infinite, because 
one can have an arbitrary relative motion of particles in a state and the 
range of relative momentum is infinite. The corresponding subspace in 
the Fock space is also unlimited. This happens despite that a regulated 
Hamiltonian cannot change a relative momentum by an arbitrarily large 
amount. The reason is that the unknown wave function ■0p(l) ■■■,n) in Eq. 
(C.6) does not limit the magnitude of relative momenta. For example, a 
square-integrable function, such as Gaussian, quickly falls off as a function 
of relative momentum of two particles but the range of relative momentum 
is infinite. 

In order to introduce a space of states with a finite range of relative 
momenta, we impose a cutoff A on the invariant mass. Namely, we define 
the space of states that are in R^pN and whose invariant mass < A. 
Ai is calculated using kinematical momenta and eigenvalues of T-Lq. The 
resulting space with this cutoff is denoted by Rat-iPN, or shortly R. We 
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introduce the projection operator on the space R and denote this projector 
also by R. Hamiltonian H\i~f projected on the space R is denoted by 

'HxNR = RT-LxnR- (C.7) 

The successive approximation to T-Lx order N in appHcation to states of 
the type \ip) with accuracy to r actions of T-Lx on with invariant mass 
cutoff A is obtained by solving equation 

Knr = [['Ho,R'HxNpR], R-HxnR] , (C.8) 

with initial conditions specified by the canonical Hamiltonian of LF QCD 
with counterterms. 

Structure of the counterterms is found by demanding that all matrix 
elements of the effective Hamiltonian Hxnr for a large ratio A/ A among 
basis states with invariant masses much smaller than A do not depend on 
the regularization applied in the canonical Hamiltonian in the limit of this 
regularization being removed |73) . 

The degree of dependence of the counterterms on N and r that are 
used in defining an approximation order N to T-Lx using r actions on var- 
ious states IV') requires studies. Specific choices of states \ip) may simplify 
identification of structure of counterterms in the canonical LF Hamiltonian 
for QCD. However, fixing all finite parts of the counterterms will require 
systematic studies of symmetries that are expected to relate finite parts of 
different counterterms, and input from phenomenology in order to fix mass 
parameters as required, in addition to the value of Aqcd- We proceed to 



showing narrowness of lixNR that satisfies Eq. (C. 
For notational brevity, Eq. ( C.9| ) is written as 



H' = [[Ho,Hjp],H] , (C.9) 

where 

H = RUxnR, (C.IO) 

Hq = R-HoR, (C.ll) 

Hi = H - Ho = R UixnR , (C.12) 

Hip = RUixNpR- (C.13) 



Evolution in A according to Eq. (C.8), or (C.9), preserves traces of 
powers of H. Trace is defined by summing diagonal matrix elements in the 
set R. The diagonal matrix elements are proportional to 5'^(0) in momentum 
space which has interpretation of volume and can be divided out. The 
argument corresponds to the conservation of total momentum. 
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Using the condition that the trace of does not depend on A, one has 
= Tr [Hi + 2HqHi + H'j)' . (C.14) 

Since Hq is fixed, 

{TrH^)' = -2TrHoHj. (C.15) 

Thus, using the basis in space R that is built from normalized eigenstates 
\m) of Hq, with eigenvalues Em, so that Ho\m) = Em\m) and Hmn = 
{m\H\n), etc., we have 



(^l-ff/mnM — ~2 Em H J mm 
mn / m 



(C.16) 



But Eq. (C.9) implies that 



H'lmn - [[Ho,Hip],Hol^^ + [[Ho,Hip],Hi]^^ (C.17) 

— {Em En) PmfiH imn 

+ [{Em - Ek)P^l + {En - Ek)P^^] HjmkHlkn • (C.18) 



The momentum is the + component of the total momentum of the 



particles that undergo interaction in the relevant matrix element {i\Hj\j). 
Momenta of the spectators of the interaction do not contribute to Pj^. The 
result needed on the right-hand side of Eq. ( C.16[ ) is 

H'lmm = '^{Em-Ek) + Pkm) HimkHlkm- (C.19) 

k 

Therefore, the sum on the right-hand side of Eq. ( |C.16 ) is 

Y^EmH'imm = Y^^Ek - Emf\HikmWP^l + P^^i) /2 , (C.20) 
m km 

and 

[Y \HlmnA = - Y^Ek - Em?\Hlkm\\P^l + P^^) < . 
\ mn / km 

(C.21) 

This result leads to the following conclusion: The sum of moduli squared of 
all matrix elements of the interaction Hamiltonian decreases with A until all 
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off-diagonal matrix elements of the interaction Hamiltonian between states 
with different free invariant masses vanish. Thus also: Matrix elements 
of the interaction Hamiltonian on the diagonal and between states of equal 
invariant masses and between states of small invariant masses as measured 
by Hq may stay constant or even increase when A decreases but only at the 
expense of still faster decrease of the off-diagonal and large invariant mass 
matrix elements and only until all off diagonal matrix elements between 
non- degenerate states are eliminated. 



Another way of writing relation (C.21) is 



\y.\'Hi^nn = -2Y,{Ml^-Mlkf\'Hikm? (C.22) 

\ mn / km 

where Mab denotes an invariant mass of particles in state a that inter- 
act in an action of T-Lj once on particles in state b. This form concludes 
our demonstration that RGPEP provides effective Hamiltonians which are 
narrow in the invariant mass of interacting particles beyond perturbative 
calculus. The invariant masses result from cancellation of spectator contri- 
butions to P~ and multiplication only by of particles that are involved 
in an interaction, i.e., not including spectators. 

Initial comparison between the non-perturbative RGPEP calculus de- 
scribed in this Appendix and perturbative RGPEP calculus described ear- 



lier, can be done by solving Eq. (C.l) using expansion in powers of a cou 



pling constant, such as g\ in the case of QCD, and observing how results of 
one way of calculating compare with the other for terms order gx and g'j^. 
Using notation introduced earlier (e.g., see [2S|, Sec. III.C, or [lOOj, Sec. 



II. B), one rewrites Eq. (C.l) as 

= -ac^ [ni]ac + Y^iPab ab + pcb cb) [Hi Ui],, , (C.23) 

b 

where letters a, 6, and c, denote configurations of particles. Using expansion 
n = Ho + 5^1 + 5^H2 + ... , (C.24) 
one obtains for the first two terms equations 

-H'lac = -ac'niac, (C.25) 

Kac = -ac^n2ac + ^{pabab + pcbcb) [UiUi]^^ . (C.26) 

b 

Using the particle size parameter s = 1/A, one obtains the solutions 

-Hlae = e-"^'^*-Hlac(0), (C.27) 
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b 



2^4 Pba ha + Pbc he r_(a62 +fec2-ac2)s4 _ 

ha? + bc^ — ac^ L 



(C.28) 



The first-order result is Gaussian in the change of invariant masses squared, 
identical to the first-order result in perturbatively defined RGPEP. The 
second-order result matches the perturbatively defined RGPEP, e.g., see 
Eq. (2.22) in Ref. jlOOj . when ac can be neglected in comparison to ah or 
be. When ac refers to a dominant Fock component and ba and be refer to a 
component with additional constituents, this condition may be satisfied very 
well if the additional constituents introduce a significant contribution to the 
invariant mass of the subsystem that counts. Such situation is encountered 
in the dynamics of heavy quarkonia [18j, where the additional constituent 
is an effective gluon while kinetic energies of quarks cannot change by much 
because quarks move slowly with respect to each other and their masses are 
large in comparison to Age* D ■ This means that the non-perturbative formu- 
lation of RGPEP can be directly applied to heavy quarkonia with arbitrary 
motion as a bound state in QCD, and implies that the non-perturbative for- 
mulation can be used to derive corrections to the harmonic oscillator force 
in Ref. [18j beyond perturbation theory to describe the spectrum of charmo- 
nium and bottomonium states that results from gluon dynamics neglecting 
light quarks. 

The procedure described here differs from truncation in powers of a 
bare coupling constant g, truncation in the number of Fock sectors, and 
combinations of both these truncations. Although we use the Fock space 
subspace R for solving Eq. (C.8), and it is useful to use an expansion 
for Tix in powers of the coupling constant g\ [16], the solutions provide in 
principle not just the matrix elements of Hamiltonian terms in a limited 
space of states but the matrix elements of Hamiltonian terms from which 
one extracts the coefficients hx of particle operators in the Hamiltonian 
terms that a priori act in the entire Fock space. Neither the Fock space nor 
perturbation theory limitations prevent us from seeking a non-perturbative 
solution for the coefficients hx. 

Some comments are in order here. One can alter the generator in Eq. 
(C.8) by introducing convergence-improving factors that are required for 
perturbative approximations [45', '36] . For bound states of low-mass hadrons 
alone, the relative motion of constituents can be described using other basis 
functions than plane waves, such as the oscillator basis. Scattering states of 
hadrons may require a more elaborate basis in coupled channels. Reasoning 
described in this Appendix justifies methods used in nuclear physics with 
only kinetic energy in SRG generator |10H 11021 11031 11041 llOSj . including 
the impact of bound states |106) . The impact can only lead to increase of 
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interaction matrix elements that are small since the combined strength of all 
interaction matrix elements must only decrease until it reaches a minimum. 
This result extends the range of applicability of SRG with a kinetic energy 
or a similar term in the generator. It suggests applications of simpler flow 
equations than Wegner's in relevant areas of physics. 



Appendix D 

Visualization of RGPEP scale dependence of hadronic structure 
Quantum fields for effective particles of size s = 1/A in RGPEP, see 



Eqs. (39), (40), and (199), are the degrees of freedom from which one can 
build effective Hamiltonian densities. An effective Hamiltonian at scale A 
describes a hadron only in terms of quarks of size s = 1/A. Therefore, 
when the scale parameter A approaches Ac ~ Aqcd, the size of quarks, 
as measured by the strong interaction range in interaction vertices |27j . 
becomes comparable with 1/Aqcd- This in turn means that the effective 
quarks become as large as the whole hadron. Such quarks must nearly 
overlap and this is how they form a white mixture. There is some imbalance 
of color on the boundary, perhaps able to couple to vr-mesons that form a 
cloud around the quarks. Visualization of these circumstances is provided 
in Figs. [3]and|4j When A is near but greater than Ac, the constituent quarks 
become smaller and part of the color-active medium inside a hadron is no 
longer fully overlapped by quarks. Therefore, when A > Ac, there must be 
additional matter density involved in filling the volume of a hadron. 

One can use an analogy with swarms of bees. For example, a baryon can 
be seen at scale Ac as built from three swarms. One swarm is made of red 
bees, one of green, and one of blue. One of the three swarms corresponds to 
one constituent quark. There is also a fourth swarm, built from bi-colored 
glue bees. But the four swarms overlap and in nearly every part of the 
baryon volume one has one red, one green and one blue bee, or glue bees. 
The regions where the three quark swarms are overlapping are locally white. 

When one changes the scale A = Ac to a nearby A > Ac as described 
in Section |5.1[ the size s of each of the three swarms becomes smaller than 
Sc = 1/Ac, but the medium around the quarks of size s remains filled with 
the glue and quark bees that can still balance to white, too. When one 
decreases A below Ac, the mass of the glue component may decrease toward 
or stabilize while the constituent quarks fully saturate the volume of a 
hadron. 

In distinction from the parton model |TT] and related scaling pictures 
\107\ [T08j . the new element of this visualization is that the constituent 
swarms are much larger in size than the distances between their centers. It 
seems to the author that this feature of RGPEP deserves a visualization 
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Figures. The RGPEP image of a meson seen at A near Ac- 




Figure 4. The RGPEP image of a nucleon seen at A near Ac- 



because none of the visualizations of hadrons known to the author has this 
feature. Namely, common visualizations represent a baryon as made of three 
constituent quarks that fly like apples in a bag, sometimes accompanied by 
some strands of glue in some kind of a cavity of unspecified nature. Images 
associated with an increase of momentum scale typically include more tiny 
objects in two ways: either as finer pieces scattered around or as interaction 
processes among quarks and gluons (such as on the cover of Ref. [109]). It 
is hard to visualize the origin of constituent quarks and CQM potentials in 
these ways. 
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In the context of RGPEP, the analogy with swarms of colored bees 
provides the following visualization for the mechanism of formation of po- 
tentials. We start with a meson built from a quark and an anti-quark. A 
quark has opposite color charge density to anti-quark. When the effective 
particles are completely overlapping like two densities of opposite charges, 
their state is locally neutral. This picture oversimplifies the non-Abelian 
picture to the Abelian one but still offers some intuition. Imagine now that 
the centers of the initially overlapping swarms move a little bit apart. It 
is well known that the Coulomb force between two spherically symmetric 
distributions of opposite charges grows linearly with the distance between 
their centers for as long as this distance is small in comparison to the in- 
dividual sizes of the distributions. Therefore, the potential energy of two 
nearly overlapping spheres is described by a harmonic oscillator potential. 
It should be stressed that the distance between the centers of the swarms 
is much smaller than the size of each swarm. In the case of baryons, one 
needs to consider three centers of three overlapping swarms being relatively 
close to each other. 

Of course, this visualization has many drawbacks: it is classical, Abelian, 
non-relativistic, and does not provide any insight concerning the difference 
between quarks and gluons. But it does model the transformation W that 
increases the number of quarks and gluons of decreasing size s when A 
increases. These smaller quarks and gluons form the medium that provides 
the expectation value of the gluon field that is interpreted in this article as 
a gluon condensate. 
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